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Abstract —This is a tutorial introduction to transient electromagnetics, 
focusing on direct time-domain techniques. We examine physical, 
mathematical, numerical and experimental aspects of time-domain 
methods, with emphasis on wire objects excited as antennas or scat- 
terers. Numerous computed examples illustrate the characteristics of 
direct time-domain procedures, especially where they may offer advan- 
tages over the more familiar frequency-domain techniques. These ad- 
vantages include greater solution efficiency for many types of problems, 
the ability to handle nonlinearities, improved physical insight and 
interpretability, availability of wide-band information from a single 
calculation, and the possibility of isolating interactions among various 
parts of an object using time-range gating. 


I. INTRODUCTION 


broadly and qualitatively defined to include all non- 
monochromatic EM phenomena, has interested scientists 
since Maxwell’s equations were formulated. However, our 
ability to obtain analytical or experimental results in this field 
is relatively recent. The two most important reasons for this 
new ability are first, the digital computer, which has made 
possible advanced computational and analytical work, and 
second, technological developments in short-pulse hardware. 
These two factors have significantly increased the availability 
of transient results for an expanding variety of EM problems. 
Besides the benefit of satisfying a natural curiosity, one 
might ask what is the value of obtaining transient EM informa- 
tion. For one thing, such information should lead to more 
effective exploitation of short-pulse hardware in practical 
applications. Examples of these could include target identifi- 
cation via short-pulse radar, and wide-band digital communica- 
tion systems. Another class of applications arises from a need 
to assess the transient EM vulnerability of electronic systems 
to such threats as nuclear electromagnetic pulse (EMP) and 
lightning. A third motivation for seeking transient EM infor- 
mation is the real possibility for achieving a more fundamental 
understanding of EM in general. As will be demonstrated 
below, many phenomena are better understood in the time 
domain in spite of the fact that equivalent information is con- 
tained in a frequency-domain representation of equivalent 
bandwidth. And, finally, where broad-band frequency-domain 
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information is sought, it can be shown that a direct time- 
domain solution will be more efficient for certain types of 
problems. 

This introductory, and necessarily sketchy, discussion of 
transient EM emphasizes engineering aspects of transients. It is 
hoped that the material will not only inform the reader about 
the current status of this developing technology, but also 
motivate him to exploit it for his own applications. 

The paper is organized as follows. First, we define our terms 
and provide an overview of the transient EM areas to be 
covered. Next, we briefly discuss the techniques available for 
realizing transient behavior in EM. Then we consider the 
physical, mathematical, numerical and experimental aspects of 
transient analysis. Finally, we illustrate features and appli- 
cations of transient EM. We demonstrate the variety of ways 
to characterize transient behavior for a given problem, as well 
as the variety of problem types that can be handled, and 
emphasize the physical interpretability of phenomena through 
transient analysis. 

The presentation is intended to be a tutorial overview rather 
than a comprehensive review. Thus the authors have limited 
the examples given to their own work, some of which were 
developed especially for this paper. A fairly extensive bibliog- 
raphy is included, however, to provide a representative sam- 
pling of the literature. 


II. DEFINITION OF TERMS AND OVERVIEW 


This particular discussion is limited to scatterers and 
radiators in linear, time-invariant media. Transient informa- 
tion for such problems can be obtained from transformed 
frequency-domain data or from a direct time-domain solution. 
We focus on the latter area, but refer to frequency-domain 
work to clarify the contrast between the two approaches. 

More specifically, we concentrate on direct time-domain 
solutions for wire objects obtained from an integral-equation 
treatment. This narrow scope allows us to treat a few topics 
in depth, rather than many subjects superficially. Although 
this presentation is fairly long, it covers only a small part of 
the developing area of transient electromagnetics. For more 
detail, the interested reader may consult various references (in 
particular, {1]-[4]) that provide an overview of the topic. 

The terms ‘frequency domain” and ‘“‘time domain” charac- 
terize the analytical or experimental procedure used to obtain 
the EM response desired. A frequency-domain procedure 
factors out the time dependence via an eivt multiplier. A 
time-domain procedure, on the other hand, treats time as an 
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explicit independent variable. The distinction can become 
blurred in pulsed-continuous-wave or swept-frequency systems, 
but in general the frequency domain uses monochromatic or 
continuous-wave excitation, and the time domain uses impul- 
sive excitation. In the discussion which follows, we differenti- 
ate in our use of the terms “transient” and “time-domain” by 
limiting the latter to situations or procedures wherein the 
data (computed or measured) is obtained directly as a function 
of time, whereas the former includes all time-dependent data 
however it has been obtained. 

When seeking transient information for a linear system, then, 
one must decide whether to employ a frequency-domain or 
time-domain approach. Generally, early work in transient 
analysis was based on a frequency-domain formulation 
because closed-form time-domain solutions are almost im- 
possible to obtain. Even a frequency-domain approach is 
difficult to solve because of the extensive computational effort 
involving many frequencies required to evaluate a problem 
and transform it to the time domain. Initial transient results 
in acoustics [5] and EM [6], [7] were for two-dimensional 
problems because of their simplicity relative to the more de- 
manding three-dimensional geometries. Rheinstein [8] later 
solved a three-dimensional transient EM problem (the conduct- 
ing and dielectric spheres) in the frequency domain, and nu- 
merous examples resulting from EMP studies have been 
developed [9], [10]. 

The first time-domain approach to EM transient problems 
was based on physical optics to obtain the approximate 
backscatter-impulse response of a flat plate and spheroid 
{11], [12]. This work was later extended to other geometries, 
such as the cone sphere. Subsequent direct time-domain work 
has concentrated on integral-equation techniques, the primary 
subject of attention below [13]-[19]. 

Some advantages to be discussed below of these direct time- 
domain solutions over frequency-domain treatments of tran- 
sient problems are: 


1) greater solution efficiency for many types of problems; 

2) more convenient handling of nonlinearities; 

3) improved physical insight and interpretability; 

4) availability of wide-band information from a single 
calculation; 

5) opportunity to isolate interactions, using time-range 
gating (e.g., pulse reflection from wire ends, bends, etc.); 

6) possibility for more directly and efficiently obtaining 
singularity expansion method poles. 


These benefits require some tradeoffs, however. Foremost is 
the generally greater complexity of a time-domain code relative 
to the equivalent frequency-domain version, with a resulting 
increase of difficulty in its development and use. The com- 
puting capability required can also be significantly larger. 


III. TECHNIQUES FOR OBTAINING TIME-DOMAIN RESULTS 


Outlined briefly below are mathematical, numerical, and 
experimental techniques for obtaining time-domain results. 
In subsequent sections, we will explore some aspects of these 
topics in greater detail, particularly the numerical area, and 
introduce before that only enough analytical background to 
support the numerical procedures that are developed. We 
mention this to emphasize our intention to concentrate on 
achieving physical understanding rather than mathematical 
depth and detail which can be found elsewhere. 
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A. Mathematical 


1) Rigorous, Closed Form: Unfortunately, few EM transient 
solutions can be expressed in closed, explicit form in terms of 
standard functions. In spite of that limitation, some success 
has been achieved for a variety of problems. Wu [20] has 
worked out a time-domain solution for a step-excited infinite 
cylindrical antenna. Felsen [21] has developed a closed- 
form solution for the diffraction of a pulsed field from an 
arbitrarily oriented dipole by a wedge. The latter problem 
belongs to a class of scattering problems for which pulse solu- 
tions take on a simpler form than do time-harmonic or 
frequency-domain solutions [22]. Franceschetti and Papas [3] 
developed a heuristic description of the general properties of 
transient radiation from elementary sources. Analytical studies 
for a sphere [23] and wire [24] have also been performed. 
But even those time-domain solutions that can be given in 
closed form may require extensive computation to obtain 
numerical results. Consequently, most time-domain solutions 
inevitably involve substantial computer processing. 

2) Asymptotic: Asymptotic techniques may involve either 
low- or high-frequency characteristics of the frequency-domain 
approach, or the corresponding late and early responses of the 
time-domain approach [12]. In either case, asymptotic 
techniques attempt to exploit what is analytically deducible 
about time or frequency behavior for the limits indicated. For 
example, we can show in the frequency domain that the 
radiated far fields produced by a finite-sized radiator or scat- 
terer must vanish as w approaches 0. Therefore, we know in 
the time domain that the time integral of the radiated far-field 
waveform must also vanish, as that integral is what gives the 
«w = 0 component in the frequency domain. 

Low-frequency asymptotic results may also be based on the 
Rayleigh law of scattering, wherein the fields go to zero with 
decreasing frequency as w?. In the high-frequency limit, we 
might employ physical optics, geometrical optics, or the 
GTD to obtain the asymptotic behavior. 

3) Hybrid: We normally associate transformed frequency- 
domain data with real frequencies; ie., the w in e“* is a real 
number. However, we can express more generally the 
frequency-domain transfer function of a given problem in 
terms of a complex frequency. 

The singularity expansion method (SEM) [4], [25] exploits 
a special feature of scattering and radiation from three- 
dimensional objects: the simple (i.e., first-order) poles their 
transfer functions may possess in the complex-frequency 
plane. If we know the locations sy and amplitudes (or residues) 
Ra of these poles, we can easily construct a transient response, 
which is simply a series of exponentially damped sinusoidal 
terms, one for each pole. Complex-frequency-domain analyti- 
cal techniques were used initially to find the poles [26], [27]. 
More recent work shows that the poles are extractable from 
time-domain data [28], [29] or from frequency-domain data 
taken on the jw axis [30], [31]. Because SEM provides a 
simple relationship between the frequency domain and time 
domain, we regard it as a hybrid technique. 

Another hybrid procedure, called the frequency-augmenta- 
tion technique [32], incorporates elements both of asymptotic 
analysis and numerical computation. It uses interpolation to 
smoothly connect between computed low frequency and 
Tesonance-region data and geometric optics, high-frequency 
data. A subsequent transform to the time-domain of the 
wide-band transfer function thus obtained produces the 
sought-for transient response. 
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B. Numerical 


EM time-domain numerical analysis naturally proceeds from 
time-dependent Maxwell’s equations. When we express those 
equations in differential form, we can solve them in terms of 
finite-difference! approximation, sampling the unknown fields 
in both space and time. Imposed spatial boundary conditions 
and temporal initial values, which are also sampled in a 
weighted sense, lead to a system of equations which can be 
solved for the sampled values over the space considered at a 
sequence of time steps. Alternatively, we can integrate 
Maxwell’s equations, using an appropriate (usually infinite- 
medium) Green’s function. Imposed spatial boundary condi- 
tions and temporal initial values now lead to an integral equa- 
tion in which the unknowns are the sources induced on the 
surfaces over which the boundary conditions are applied. A 
sampling of these unknowns in space and time and a weighted 
sampling of the boundary values again lead to a system of 
equations for the space-time sampled values of the unknown, 
which is solvable as an initial-value problem. 

There are some important differences between these two 
approaches. First, in the differential-equation formulation, 
the unknowns are sampled at all points within and on the 
boundary(ies) of the solution space. In the integral-equation 
formulation, on the other hand, the sampling is done only over 
the boundary on which the boundary conditions are applied. 
Thus, the integral-equation method can result in substantially 
fewer unknowns. 

Second, the integral-equation treatment requires a space- 
time-dependent Green’s function. This function allows surface 
sampling to replace the volume sampling of the differential- 
equation formulation. For lossless, nondispersive, homoge- 
neous, linear, isotropic media, such a Green’s function is 
readily obtained, but otherwise substantial complications 
occur that can require a volume, rather than a surface, integra- 
tion. Consequently, for these more general media, an integral- 
equation approach may not be the most suitable: Nevertheless, 
most EM time-domain analysis to date uses the integral- 
equation approach, which is the only one we discuss in detail 
in this paper. 


C. Experimental 


Transient-response data are measured primarily through 
direct time-domain procedures. Developments in short-pulse 
technology enable us to generate and measure high-amplitude 
(103 V), fast-rise-time (~300 ps) pulses [33], [34]. By using 
these pulses to excite a test object such as a scatterer or 
antenna, and using a sampling oscilloscope or transient digitizer 
to measure induced currents and scattered fields, instantaneous- 
measurement bandwidths of 10:1 and more are possible. 
Such data can validate time-domain calculations directly; 
transformed to the frequency domain, they can meet a variety 
of needs. 


IV. PHYSICAL ASPECTS OF TRANSIENT ANALYSIS 


Perhaps the single most useful aspect of transient analysis is 
the opportunity it provides for more clearly depicting and 
interpreting the physical behavior of electromagnetic fields. 
For example, a short pulse propagating on an open-ended wire 


1 We include any method which approximates differentials by discrete 
samples in the finite-difference category, e.g., the finite-element tech- 
nique, unimoment method, etc. 
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clearly demonstrates the effects of radiation damping, disper- 
sion, and reflection from an impedance discontinuity (the wire 
end or a discrete impedance load). Because such effects are 
harder to interpret as a function of frequency, they must 
usually be indirectly inferred in the frequency domain. Of 
course, by transforming frequency-domain results, we can 
derive the transient behavior which illustrates such phenomena, 
but direct development in the time domain is generally more 
tractable, and the concepts are easier to visualize. 

In this section, we discuss fundamental physical aspects of 
transient radiation and scattering from a heuristic viewpoint. 
We first present typical results of time-domain computations. 
Then we examine in some detail the phenomena thus illus- 
trated. Next we consider radiation fields from the viewpoint 
that they are a manifestation of charge acceleration. 


A. Some Numerical Resuits 


The calculations presented here precede a detailed account 
of their mathematical and numerical aspects in the next two 
sections. They are included now to provide a mental image of 
some physical aspects involved. 

Figs. 1 and 3, respectively, show current and charge distribu- 
tions on a straight wire of length ZL = 1 m excited by a voltage 
source at its center as an antenna, and excited by a tangential 
electric field of a normally incident plane wave as a scatterer. 
In both cases, the time variation of the exciting source is 
Gaussian, i.e., exp (-g?¢7), and space distributions of current 
and charge are shown for several instants of time. Figs. 2 and 4 
show in part (a) the time variation of the current at the wire’s 
center and in part (b) the resulting far fields in the broadside 
direction corresponding to the antenna and scattering cases of 
Figs. 1 and 3. 

In Fig. 1 the current for the antenna divides into two pulses, 
approximately Gaussian in form and with a small oscillatory 
undershoot, which propagate outward from the source, ac- 
companied by oppositely signed charge pulses. The pulses 
(current and charge) decrease in amplitude and spread out as 
they propagate from the source. Nearing the end of the wire, 
the charge pulses increase in amplitude, while the current 
pulses decrease, falling to zero amplitude at the end. After the 
reflection, the amplitude of both pulses diminishes, and the 
current reverses sign because the charge flow reverses direction. 
Initially the radiated field closely resembles the Gaussian 
shape of the exciting pulse, with which it coincides in retarded 
time (allowing for propagation time) because the antenna 
behaves as though it is infinitely long prior to end reflection, 
but then exhibits a slight negative undershoot. A large negative 
pulse in the radiation field coincides with the end reflection 
of the current-charge pulse. 

The current for the scatterer distinctly differs from that for 
the antenna. It is uniform over the entire wire except near the 
ends, where it falls to zero, and where the charge is first con- 
centrated. The current amplitude in the central region of the 
wire slowly decreases with time, while near the ends it reverses 
sign as the boundaries of the pulse collapse inward due to the 
end reflection. The initial return of the scattered (or re- 
radiated) field is similar in shape to the Gaussian exciting pulse 
and coincides in retarded time with the current buildup. This 
part of the scattered field is sometimes referred to as a specular 
flash. A sign reversal of the scattered field closely follows, 
with a peak value less than half that of the first maximum, 
beyond which the field decays with time. This part of the 
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Fig. 1. Current and charge distributions at several instants of time for a 
straight-wire antenna excited by a Gaussian pulse. The antenna me 
numerical parameters are: length L = 1 m, radius a = 6.738 X 10°? 
(which gives 2 = 2 In Le = 10); Gaussian-pulse parameter g = 5.556 x 
107'° where the pulse time variation is exp [-g’f’ ]; space-segment 
length A = L/60; and time interval 6 = A/c. After its initial excitation, 
the charge-current pulse splits into two oppositely propagating pulses 
of oppositely signed charge, resulting in current pulses having the 
same sense. A slight decrease in the amplitude of the pulses can be 
observed as they progress down the wire, and a more noticeable 
decrease in amplitude occurs upon end reflection. 


Source current — mA 


Radiated fields times distance — V 


if) 10 20 30 
Time — ns 


(b) 


Fig. 2. (a) Feedpoint current as a function of time. (b) Broadside- 
radiated electric field as a function of time for the straight-wire antenna 
of Fig. 1. Six hundred time steps are shown. The peaks in the 
radiated field appear to coincide with the initial excitation and end 
reflection (in retarded time). 
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(a) () 7 (c) 


(g) (h) (9) 


Fig. 3. Analog of Fig. 1 for the scattering case. The excitation is now 
a Gaussian-pulse plane wave incident from broadside. A uniform cur- 
rent is initially excited along the wire and collapses inward due to end 
reflection. 


Load current — mA 


Scattered fields times distance — V 


0 10 20 30 
Time — ns 


(b) 


Fig. 4. (a) Transient current at the center of the long-wire scatterer. 
(b) Scattered field as a function of time for the straight-wire scatterer 
of Fig. 2. The scattered field exhibits a specular-flash return as the 
current is first excited, then a longer lasting but slowly decaying 
portion of opposite sign. 
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scattered field has a time variation similar to that of the decay- 
ing current. 


B. Interpretation of Numerical Results 


The relationship observed between the antenna current and 
radiation field clearly shows that the onset of current (and 
charge) flow is responsible for the first portion of the radiated 
waveform. The subsequent sequence of oppositely signed 
pulses in the field are as clearly due to the reflection of the 
charge from the wire ends (see Fig. 24(a) for further evidence 
of this). These results suggest that a radiated field is produced 
both when charge is accelerated (during turn-on of the exciting 
field) and when charge is decelerated (during reflection from 
the wire ends). Note that the sense of the radiated field 
evidently depends upon the direction of the charge accelera- 
tion; the deceleration here corresponds to negative acceleration 
with respect to the original charge motion, and produces a 
field of opposite sign at the first end reflection. The inter- 
mediate part of the radiated field between these first two 
pulses is not so easily accounted for. It appears to correlate 
with the spreading out of the charge pulse and with a corre- 
sponding decrease in its amplitude (dispersion) as it propagates 
down the wire, a phenomenon which can be detected if the 
results of Fig. 1 are carefully observed. Thus it could be 
ascribed to a shedding of energy during propagation. 

The initial pulse in the scattered field is also clearly due to 
the initial charge acceleration associated with arrival of the 
exciting (or incident) field. The nearly immediate sign change 
of the scattered field is evidently due to charge deceleration at 
the wire ends upon onset of reflection. In this case the result- 
ing opposite-sign radiation persists, although its value decreases 
with time, because current excited along the entire wire must 
eventually propagate to either end (recall that the antenna 
current divided into two oppositely propagating pulses), where 
it is reflected. The decreasing amplitude of the scattered field 
in this portion of its waveform must be due primarily to the 
decrease in time of the charge reaching the ends, as can be in- 
ferred from the diminishing current amplitude in the central 
area of the wire observable in Fig. 3. This decrease in current 
amplitude provides further, although indirect, evidence that 
the charge radiates as it propagates along the wire. 

In both the antenna and scattering cases, then, we find 
evidence that accelerated or decelerated charge causes radia- 
tion.? This radiation occurs whether the charge is driven by 
an accelerating, exciting field which provides an “external” 
force, or by a decelerating, induced field due to the charge 
itself which provides an ‘‘internal” or self force. The external 
force does work in setting the charge in motion and produces 
the associated radiation fields. The internal force, due to 
charge buildup at the wire’s end, receives energy from the 
charge as it slows to zero velocity. As the charge accelerates in 
the opposite direction, not all the energy returns to it. The 
difference is lost in the form of radiation. 

An additional loss seems related to dispersion of a charge 
pulse propagating on a wire, a result directly observable in the 
antenna case and inferrable in the scattering case. The loss 
mechanism in this case might be the charge deceleration which 
the pulse dispersion implies. We might conclude that this be- 
havior arises because the trailing edge of the pulse propagates 
more slowly than the leading edge, whose speed is near that of 


2In the following, we use the term acceleration to cover both cases. 
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a light wave in the medium. Propagation of the charge pulse 
in this situation might be likened to the behavior of electron 
bunches in a klystron. It should be understood, however, that 
charge flow on a conductor differs from charge flow in space. 

In the latter case it can be a single charge in motion that is 
the source of a pulse of EM fields propagating through space. 
In the former case, a thin wire for example, a charge-current 
pulse like that of Fig. 1 is not comprised of the same charges 
all along the wire. Instead, a local charge separation is caused 
in response to the fields other charges produced at earlier 
times, giving the effect of a moving charge pulse. It is possible, 
therefore, to consider the movement of charge along a conduc- 
tor as equivalent in some respects to charge movement in free 
space. 


C. Some Mathematical Relationships 


It is indeed a fact that charge acceleration causes electro- 
magnetic radiation. This fundamental principle provides a 
basis for understanding all radiation processes, both transient 
and steady state. We emphasize the role of charge here, since 
while charge and current may both appear as source terms in 
Maxwell’s equations, charge can be viewed as the more funda- 
mental quantity. Current, on the other hand, is a derived 
quantity resulting from the charge motion. 

We can study accelerated-charge fields in two ways: the 
microscopic approach, which concerns fields of individual 
charges, and the macroscopic approach which concerns average 
fields over the charge distributions [35]. We primarily con- 
sider the latter approach here. We can distinguish between the 
two on the basis of whether the observation times and dis- 
tances are smaller (microscopic) or larger (macroscopic) than 
the characteristic times and distances associated with the 
sources [3] (e.g., the relaxation time in a metal or its skin 
depth). The microscopic approach involves equations of 
motion of individual electrons, and can include mass, relativis- 
tic, and quantum-mechanical effects. In the macroscopic 
approach, we are more interested in solving Maxwell’s equa- 
tions in a field description than in the physical details of their 
sources. 

A macroscopic viewpoint of radiation is provided by the 
integral form of Maxwell’s equations from which the total 
fields due to an electric current distribution J and charge 
distribution p over a perfectly conducting surface S can be 
written 


WF',7) 1 8 = R 
- 24.5 @n}x oy ds’ (1a) 
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The source and observation point locations are at 7’ and 7, 
respectively, ¢ is the observation time and 7, the retarded time, 
is t- R/e with R=7- 7 and R=|R|. Also, uo and ég are the 
permeability and permittivity of free space, respectively, and c 
is the speed of light. In the far field, where |7] =r > %° and 
Rr these equations simplify to 
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Eraa (7 t) = No H(, t) X = (b’) 
where to, the wave impedance, is given by Vo /Eo - 

Since by definition J is the flow of charge, and charge is 
conserved, any change in J with time is equivalent to a change 
in velocity of the charge which carries the current. Equation 
(la), therefore, states that radiation is due to accelerated (or 
decelerated) charge, in agreement with the preceding discussion. 

A microscopic viewpoint of radiation is provided by the far- 
radiated fields of a point charge q moving with (nonrelativistic) 
velocity v << c¢ and acceleration v which are [36] 


we x 4 furyXR  RXRX [RX i))I 

a, o£ R3 - cR* a 
_ qa {R RXRXiG@) 

s.o= 45 a al (2b) 


When v is zero, H and E have only 1/R? static components, 
an elementary result, and Exad changes sign when v reverses 
direction. The rate of energy radiation is also given by 

dw 2;=)2 

ati Ne (2c) 

dt 67Egc 
The total amount of energy radiation during a given velocity 
change is proportional to the time integral of |v|?, so the 
higher the acceleration, the greater the radiated energy. 

In our time-domain solution, we can also develop alternate 
expressions which provide a macroscopic description of the 
radiation process in terms of the current and charge distribu- 
tions on the object. For the particular case of a wire object 
having contour C(7), the net energy collected by the object 
from the exciting electric field FE... up to time ¢ is given by 
(see Appendix) 


t 
Wo = W(t) | collected -{ i Eexc(w’, t') 
-0 “C(7) 


-PI(w',t')dw' dt. (3a) 


Similarly the total energy dissipated due to resistive loss R(w’) 
can be expressed up to time ¢ as 


t 
Wp = W(t) |pissipated = f Pw',t')R(w') dw’ dt’. 
- YC(r) 


(3b) 


Finally, a measure of the total energy stored in the fields near 
the wire due to its current and charge density at time f, is 
provided by 

1 2(w' t 
Ww = W(t)lwire * — [uote t+ Pon) dw’ 
c@) £0 


=W, X 10°? + Wg X9X 10° (3c) 


where W; and Wo are the current and charge contributions, 
respectively, with J(w,t) and Q(w, ft) the current and charge 
space-time distributions, and # the unit tangent vector to 
C(F) at 7. 

Results for W; and Wo as defined by (3c) are presented for 
the radiation and scattering cases just considered in Figs. 5 
and 6, respectively. Referring to Fig. 5 for the Gaussian-pulse 
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Fig. 5. Long-wire antenna results for the case already considered in Fig. 
1. Wy and Wo are derived from the integral over the antenna of J? and 
Q’, respectively, and thus represent a measure of the energies stored 
on the wire in the current and charge. End reflection of the current- 
charge pulse (see Fig. 1) results in a substantial reduction of these 
quantities. 
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400 


Fig. 6. Current and charge energy Wy; and Wo stored on a long-wire 
scatterer for the case already considered in Fig. 3. 
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excited antenna, we observe several interesting features. In the 
smoother parts of these curves (away from the sharper peaks 
and nulls), we find that the ratio of W; to Wg is ~c? , showing 
the current and charge contributions to the total energy are of 
equal value, according to (3c). Also, in these smoother 
regions, a decrease in time of both W; and Wg occurs due to 
radiation damping of the propagating current-charge pulses, a 
result already inferred from the current plots of Fig. 1. At 
times of end reflection, Wy; tends toward zero because the 
current vanishes at the end while Wg is sharply peaked, show- 
ing the energy storage is then primarily in the charge. Con- 
versely, when the charge pulses meet at the wire’s center, Wo 
becomes zero due to charge cancellation and W; is peaked. 

For the scattering case of Fig. 6, the results are not quite so 
readily interpretable. One thing is clear, however, and that is 
the periodic interchange which takes place between W; and 
Wg. Observe that the maxima and minima of one are co- 
incident with the minima and maxima of the other, and that 
their sum decays monotonically in time. 


D. The Role of Current and Charge 


In the antenna case, the excitation is a tangential electric 
field applied at the center of a small region of length A (of 
which there are N where AN=Z) on the antenna’s surface. 
This region for antennas is commonly referred to as the ‘‘gap,”’ 
(or source region) the term we employ below, although it 
should not be interpreted to be a physical break in the 
antenna’s surface but instead only denotes where an antenna 
is excited mathematically. For simplicity, assume the field is a 
unit-amplitude pulse of time duration 6, where A=cé, ap- 
plied at time t=0. Our model is thus discretized in space 
(segment of length A) and time (intervals of duration 6) and is 
conceptually identical to the numerical solution developed 
below. 

Because the assumed perfect conductivity requires the total 
tangential electric field to be zero, initial application of the 
exciting electric field instantaneously induces a charge separa- 
tion in the gap, which results in charge acceleration and initia- 
tion of an associated current. The tangential electric field of 
this separated charge cancels the applied field. Due to the 
finite propagation velocity of the induced fields resulting from 
this charge motion and separation, there is no effect produced 
outside a region of radius A centered in the gap during the 
time 6 the excited field is applied. 

Upon removal of the exciting field after a time 6, the charge 
separation in the gap must vanish to make the total field there 
zero. Simultaneously, the current and charge disturbance 
produced in the gap during 0 <r <6 begins to affect a region 
A wide on either side of the gap during the time 6 <¢ < 28. 
This effect takes the form of a reradiated or “scattered” 
tangential electric field having the same directional sense as 
that of the original charge separation in the gap. Conse- 
quently, a compensating net charge density must occur in each 
of these regions to cancel this scattered field. The continuity 
equation shows that this is consistent with the movement of 
the charge from the gap segment to the adjacent segments (as 
shown in Fig. 1). Charge neutrality is thus restored in the gap 
and current flow there is zero. A current flow of half the 
original gap current is set up in these two adjacent regions. 
This current has half the original value because it involves the 
movement of only half the original charge (only the plus or 
minus charge, not both as in the gap). The current has the 
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Fig. 7. Radiation can originate from several places on a wire object. 
Here they are depicted conceptually and include: (a) the source 
region; (b) an impedance load; (c) a change in radius; (d) sharp bend 
(or a junction); (e) a smooth curve; and (f) an open end. Also a cause 
of radiation is the dispersion which occurs on a straight section of 
wire, which can be seen in Fig. 1 in the decreasing amplitudes of the 
current and charge pulses as they propagate along the wire. 


same sense on either side of the gap because it involves the 
oppositely directed flow of oppositely signed charges. 

After a series of repeated steps like that above, these pulses 
of current and charge reach the ends of the wire. The charges 
stop because the conduction current must vanish there. In 
stopping, the charges are subject to an oppositely directed ac- 
celeration and so produce a scattered field in a sense opposite 
to that of the original scattered field. A reflected charge and 
current pulse are thus produced, and the process described 
above continues, but with some loss due to radiation. From 
(2c), the amount of radiated energy is proportional to lvl?: 
the shorter the time interval over which a given amount of 
charge reflects, the more the radiation loss. Thus the narrower 
the incident pulse, the more efficient the radiation production 
upon end reflection. 

We can analyze the scattering problem in the same way. Two 
factors which influence field and current behavior and which 
help us to understand the overall problem include: 1) causality, 
which limits the space-time region over which interactions can 
occur, and 2)the current-continuity equation, which con- 
strains current and charge behavior. Energy is supplied to the 
object while the exciting field is present. All this energy may 
be regarded as being radiated, but only a portion actually 
produces a space wave (far-field radiation). The rest remains 
associated with the structure along which it propagates. At 
each point where energy flow is redirected by the object, 
reradiation occurs. As in the source region, part of the energy 
leaves in a space wave and part remains attached to the struc- 
ture. Where energy-flow redirection does not occur, as along a 
straight wire, there occurs no radiation (unless there is dis- 
persion), but only a continued propagation of fields previously 
produced. 


E. The Field Behavior 


A perfectly conducting object can have an electric field at its 
surface in only the normal direction. The total tangential field 
is zero by definition. Thus, without a normal field, there is no 
field at all. Consequently, for radiation to occur, there must 
be: 1)a nonzero charge density to produce a normal electric 
field, or 2) an exciting electric field that is cancelled by an 
oppositely directed induced field. The latter cause of radiation 
is charge acceleration, as discussed above. The former is 
similarly associated with charge deceleration. The radiated 
electric fields may be closed (i.e., form continuous loops) or be 
open as a result of terminating on the object. They can close 
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Fig. 8. During the time that an impulsive field excites a perfectly conducting object, there occurs a net flow of energy from the incident field 
onto the surface because  - Eex¢ # 0, part (a). Because  - Etotay = 0 due to its infinite conductivity, there is an equal and opposite flow en- 
ergy outward in the buildup induced or scattered field. This field is accompanied by buildup of an induced current and charge distribution 
which provides a measure of the energy stored in this field. After Eexc becomes zero, part (b), this stored energy is eventually radiated. 
This happens even though there is zero net energy flow at its surface. We suggest as a model that the stored energy “reflects” from the ob- 


ject’s surface. 


Some of this reflected energy is radiated or lost to the far field; for an antenna, this takes place most efficiently when the 


current charge pulses change direction (e.g., reflect from the open ends of a straight wire). That which remains is gradually lost with each 


subsequent reflection. 


only as a result of charge cancellation, caused, for example, by 
the meeting of oppositely signed charges moving in opposite 
directions on a wire. The radiation fields of finite-sized, 
perfectly conducting objects produce radiation fields with 
closed lines at distances greater than the maximum object 
dimension. 

Near fields, i.e., those field components that do not carry 
energy away from the object, include electric fields that 
terminate on the object due to its net charge density and 
associated magnetic fields that result from charge motion. 
These near fields contain energy which is essentially bound to 
the object’s surface, and which acts as a guided-wave structure 
for near-field energy propagation. Charge deceleration con- 
verts near-field energy to radiating energy. It is useful to think 
of energy due to charge separation as potential energy, and 
that due to charge motion as kinetic energy (see (3e)). Both 
the fields and charge share the stored energy. Further, the 
charge terminates the fields on the conducting surface through 
the boundary conditions that Egy and Hyorm are zero. 


F, Summary 


The following processes are primarily responsible for electro- 
magnetic radiation from wires: 


1) charge acceleration due to an exciting electric field; 

2) charge deceleration due to reflection from impedance 
discontinuities such as discrete loads and open ends and 
abrupt bends; 

3) charge deceleration due to smooth curves; 

4) charge deceleration due to dispersion. 


and are illustrated qualitatively in Fig. 7. 

We may summarize the above discussion by invoking the 
following model for transient radiation from a perfect conduc- 
tor. It is based simply on the fact that, due to the vanishing of 
the total tangential electric field on the conductor’s surface, 


the net energy flow (Poynting’s vector) across it must always 
be identically zero. During the time the exciting field (source) 
is present, the outward flow of energy due to the induced cur- 
rent and charge are thus exactly balanced by the source- 
supplied energy. After the exciting field becomes zero, there- 
fore, there can be nothing other than zero energy flow normal 
<Wy) to the conductor’s surface. But since the conductor 
must continue to radiate until all the energy is gone, we sug- 
gest that Wa has two counter-flowing components. One is an 
inward-flowing component which contains energy collapsing 
back onto the conductor from the fields near it. The other 
is an outward-flowing component which contains energy that 
propagates away from the conductor and is then lost. These 
processes are depicted conceptually in Fig. 8. 

Both components must remain everywhere equal as they de- 
cay monotonically with time until eventual charge neutrality is 
restored. Although the total near-field energy is not easily 
derivable from the current and charge, these induced sources 
provide a convenient and useful indication of the stored 
energy, and thus also of the radiation process, as we have seen 
above. 


V. MATHEMATICAL ASPECTS OF TRANSIENT ANALYSIS 


In this section we consider mathematical aspects of direct 
and transform techniques for obtaining transient results. We 
continue with a discussion of hybrid techniques that possess 
characteristics of both frequency- and time-domain analysis, 
and conclude by briefly discussing the limitations of these 
various approaches. 


A. Direct Techniques for Obtaining Transient Results 


From Maxwell’s equations for a perfectly conducting closed 
surface we can derive two time-dependent integral equations 
based upon the electric and magnetic field, respectively, as the 
forcing functions [2]. 
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Most frequency-domain antenna analysis is based upon the 
thin-wire approximation. This involves replacing a two- 
dimensional surface integration with a one-dimensional line 
integration while approximating the surface current as an 
axially directed filament, and most often uses the electric-field 
type integral equation. This same approximation is also 
useful in the time domain, and leads to 


ww’ 9 ' 
® + Bexc(7, t) = rae = 37 137) 


Raow',2| dw’, FECta (4a) 


where 


, i a , t t 
q(Ww ,m=- — Iw ,t)dt 
co OW 
% and ®’ are tangent vectors to the wire at 7 and 7’, respec- 
tively, C is the wire contour, a is the wire radius, andr EC +a 
denotes that the field is to be evaluated on the wire surface. 
Also, I(w, t) is the wire current at point w and time t, R= 
|R|=|7-7'|,and7=t- R/c is the retarded time. Other time- 
domain integral equations specialized to wire geometries can 
also be developed. Equation (4a) provides the basis for most 
of the subsequent discussion in this paper. The far-field ex- 
pression which corresponds to (4a) can be written 


2 ry! T) +e% 2100] dw' 


Fir ya - He ’ 
EG, 0) C E ar , R dw 


(4b) 


B. Transform Techniques for Obtaining Transient Results 


Frequency-domain analysis also proceeds from Maxwell’s 
equations expressed in differential or integral form, but for an 
assumed e!? time variation. Imposed conditions on the field 
behavior over spatial boundaries complete the analytic formu- 
lation of the problem. A frequency-domain counterpart of 
(4a) for a wire geometry is given by 


0 BaetF)= 2 | I(w') 
4m Jo 


[o-0+ 20-00: ]eaw FECta, 


where 


y=exp [-ikR]/R (5) 


and 


K= WV Uo€o- 


C. Hybrid Techniques for Obtaining Transient Results 

Hybrid techniques, with characteristics of both frequency- 
and time-domain analysis, include the frequency-augmentation 
technique and the singularity-expansion method discussed 
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below. Other possibilities might also be identified, such as 
various asymptotic techniques mentioned above. 

1) Frequency-Augmentation Technique: Bennett et al. [32] 
developed a technique to combine low- and high-frequency 
scattering information for an object to obtain its impulse re- 
sponse. The basic idea of this frequency-augmentation tech- 
nique is to represent the impulse response of the object in the 
frequency domain H(w) as 


A(w) = F4(w) + H4(w). 


The first term, the augmentation function, approaches H(w) as 
« approaches °, and thus represents the high-frequency in- 
formation referred to above. The augmentation function 
might come from either time- or frequency-domain analysis, 
and could be based upon geometrical optics, physical optics, 
or other considerations. The second term, the augmented- 
frequency response, thus vanishes as w approaches ©, and is 
derived from low-frequency information as follows. 

Consider, for example, that, using either the time- or fre- 
quency-domain approach, the actual low-frequency response 
function of the object is known up to some maximum fre- 
quency W,, and is denoted by A(w). Then define 


H4(w) = A(w) - Fa(w) 
an ote that 
H4(w) = H4(w), 


But since #,4(w) is known to approach zero with increasing 
frequency, we can introduce some reasonable asymptotic 
function to represent this behavior and match it to Hi,(w) in 
the region below w,, to obtain an approximation for H4(w) 
for all w. If the function so chosen is accurate enough, then 
we obtain a combined numerical-analytical representation for 
H,(w). Finally, a Fourier transform to the time domain yields 
the desired impulse response, h(t). The augmentation function 
contains the singular, high-frequency part of the response, and 
can be analytically transformed to the time domain, thus avoid- 
ing the problem of trying to numerically transform a singular 
function. Only H(w) needs to be transformed using a numeri- 
cal FFT. 

2) Singularity-Expansion Method: We have already men- 
tioned the basis of the singularity-expansion method (SEM). 
Here we very briefly discuss some of its basic characteristics. 
SEM stems from the observation that the late-time EM response 
(i.e., when the exciting field no longer is present) of an object 
has the form 


for WS Wy. 


N 
f(t) = D° Rg exp (Sq?). 


a=1 


This time function has a similarly simple complex-frequency 
representation, 


F(s)= y 


Note that f(t) and F(s) may represent any measure of the 
object’s response, including induced current, near field, and far 
field. 

In transient solutions, then, the original problem is replaced 
by the new requirement to find the Rg and sg. This require- 
ment can lead in turn to at least two distinctly different ap- 


= - e’ 
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proaches, depending on whether time-domain or frequency- 
domain data are available for this purpose. A procedure based 
upon Prony’s technique has been successfully used to obtain 
poles both from time-domain data and frequency-domain data 
{26}, [30]. 

It is extremely important to note that the s, depend only on 
the object geometry. Thus once determined, they characterize 
the object for any excitation and can be useful for various 
purposes, including radar target classification [30}, [31]. The 
Rg, on the other hand, depend on both the object geometry 
and its excitation. 


D. Limitations 


Our ability to obtain transient results is limited by the same 
factors that apply to the frequency domain. Analytical solu- 
tions are available for relatively few problems, and those that 
have been developed generally require extensive computations. 
Numerical techniques are limited primarily by what is com- 
putationally feasible. From a mathematical viewpoint, then, 
our ability to formulate transient problems greatly exceeds our 
ability to subsequently cast them in a numerically tractable 
format. However, alternate formulations may exist which 
could greatly reduce the numerical effort required for solving 
a given problem. 

Generally, the numerical approach is most effective for fre- 
quencies up to and including the resonance region, for the 
relatively simpler geometries. Any increase in complexity of 
object geometry or environment can greatly reduce or even 
eliminate the possibility of using a numerical approach. This 
is one area where further formulation efforts might be directed. 


VI. NUMERICAL ASPECTS OF TRANSIENT ANALYSIS 


The frequency domain and time domain share a common 
numerical solution procedure which is the moment method. 
Proceeding from an analytical statement of the problem (an 
integral equation, for example) to a sampling of the unknown 
sources and the known boundary values leads to a linear sys- 
tem of equations, which can be solved numerically using matrix 
techniques. Sampling of the unknown sources and known 
boundary values in the frequency domain is done only as a 
function of space, whereas in the time domain this sampling 
is done as a function of space and time. In the following dis- 
cussion we present the numerical ingredients of the time- 
domain approach, and a comparison of the computational 
characteristics of both the time and frequency-domain pro- 
cedures for obtaining transient information. 


A. Direct Techniques 


For clarity, we illustrate the direct time-domain solution 
procedure for a simpler, or prototype, time-domain integral 
equation in place of the actual (4a). We thus consider 


h 
ee. { f(x’, t')K(x,x') dx’, -hXh (6) 
-h 


where 


t(x,x', N=t- |x- x’ |/e 


as the equation to be solved, where g is specified and f is to 
be found. Proceeding on an intuitive basis, we might decide 
to approximate (6) with a discrete sequence of samples for 
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f as a function of both x and «. If we further choose some 
reasonable variation between these discretely sampled points 
(i.e., select an interpolation function) and also specify how the 
right and left sides of this sampled equation are related, then 
we can reduce (6) to a linear system in which the samples for 
f are the unknowns. 

For practical purposes this procedure constitutes the mo- 
ment method. The interpolation function mentioned, in the 
moment-method context, is called a basis function, and the 
telationship between the two sides of the reduced equation 
depends on a weighting or testing function. Let us now use a 
space-time pulse approximation for f, i.e., 


N i 
fe’, t= > SO UpyAr;' 


i's1 j'=1 
where 
xj - A/2<x' <x, + A/2 


U;';' = 1, if ; 
: oe ee one 


and is zero otherwise, and x,’ = Ai’, t= 57’. Let us also point- 
match the integral equation at the space-time sample locations 
x;, tj. Then we obtain 


8ij = 8 (x;, tj) 


h N j 
= DD Uy Agj'K(e;, x") ax! 
-h i’=1 j'=1 
N (i'+ 1/2yA j 
= >. f > Uy 'Ar;'KGA,x')dx' (7) 
fei 4(i'-12)A j’=1 
where i= 1,:-+,N andj=1,°++,Nr. But since t’=t- |x - 
x’ |/c, and upon using A = c6, equation (7) simplifies to 
N (i'+1f2)a 
Si; = >; Aj! j-\i-i'| J KA, x') dx' (8) 
i’=1 (i'-1/2)A 


which associates with each A;‘;', a spatial integral (or interac- 
tion coefficient) over segment i’ at retarded time j - |i - i’|. 
It is helpful to rewrite (8). Let us first denote the interaction 


coefficients by 
(i'+1/,2)A 
Zii' -{ K(iA, x’) dx’ 


('-1/2)4 
so that 


N 
Bi= 2) Zu Ar',j' 
i'=1 
= LZyAy + Z;,4-1Aij-1,j-1 
+ Zii-2Ai-2,j-2 °° + 271 Ay j-i91 
+ Zp i+ 1 Aisi,j-1 + 254424 i+2,j-2 


(9) 


We observe that the first term involves a sample at the time 
step 7, while all other samples are forj- 1,j7- 2,°°-,and thus 


+:° *+ZipnANn,j-Neoi- 
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are from earlier times. Upon solving (9) for A ijs We then find 


N 
By DL Za Ai j-ti-i' 
- i’=1 
Ay= Tee pee, ee 
i“ 


f=1,°°, : 
,i #i (10) 
fRle Np 


where the summation excludes the term i’ =i. 

Equation (10) should forcefully demonstrate that at time 
step j the integral (6) can be solved as an initial-value problem 
by time-stepping, for if all values of A;, are known for k <j - 1, 
then Aj; is completely specified by them and the present value 
of the forcing function g;;. This equation shows the explicit 
effect of causality and the finite velocity at which EM fields 
propagate, as an increasing time delay between a source at 
i’ and its influence being observed at i (through the time index 
j-|i-i'[). In principle this factor permits solving (6) without 
matrix inversion, as further emphasized by the Z;;, coefficient 
matrix in (10). 

We have of course simplified this discussion of the numerical 
solution by choosing an apparently simple integral equation 
and using pulse-basis and delta-weight functions. The actual 
computer code [37] used to generate the results presented 
here uses a nine-term polynomial basis, up to and including 
quadratic space-time variation, and delta-function weights. 
Nevertheless, because the integral equation used here repre- 
sents that which would apply to a straight wire, it is realistic. It 
makes retarded time depend very simply on source-observation 
distance, i.e., 


t=t-|x-x'I/c. 


Therefore, as we integrate the source space (x’), the space-time 
path described by x’ and t’ isa straight line, as illustrated in 
Fig. 9. This integration path furthermore passes diagonally 
through the center of each space-time sample “patch” (or Aj), 
because the observation points are also located at the patch 
centers and A=c6 was used everywhere. This path means that 
a single Aj; is associated with each Z;;'. In general, the inte- 
gration path can pass through two Aj,’s for integration over 
one space segment when using a pulse basis, so that a given 
Zjj, decomposes into two parts, one multiplying each respec- 
tive Aj; value. 

For example, if the space observation points are not collinear 
(i.e., on the same line defined by the current sample), then 
t' and x’ are not linearly, but instead 


t=t-Va-xP+iy-y' YP +-2' Ye 
=t- Vp? +(s-5')/c 


where s and p are cylindrical coordinates of the observation 
point relative to the line along which the source current flows. 
Fig. 9 also shows other representative space-time integration 
paths. 

The g term in (6) (or the g,; terms in (10)) represents the 
actual excitation (or its sampled values) which excites the 
response f (or the sample values A;;). Since we are dealing 
with an integral equation based on the electric field and spe- 
cialized to a wire, g represents the tangential electric field 
distribution along the contour C of the wire. For a scattering 
calculation, and assuming a point-matching solution is em- 
ployed (which is the case for all the results subsequently pre- 
sented), gj; then is simply the specified value of the incident 
electric field tangent to C at a sequence of space points r; and 
time steps t;, For an antenna calculation, the exciting field 


PROCEEDINGS OF THE IEEE, VOL. 68, NO. 11, NOVEMBER 1980 


i+] 


Observation point 
o fj) 


j e 


5s j-2 e ° . 
& Source space-time 
E patch i-4, j-4 
= 4j- 3 e e 
-! 0 i 
Alc 
ind ghee 2 schol - 4 


Integration paths in space 
and retarded time 


Space index 


Fig. 9. Space-time integration paths in the source-coordinate (s’, ¢’) 
space for an observation point located at time ¢t and cylindrical coor- 
dinates p and s relative to the center axis of the current filament (seg- 
ment) being integrated. The retarded or source time f’ is given by 
t=t- Jp? +(6-s'/ce. Note: s' = i'A;s =i; ft = j'8; t = 8. The 
uniformly spaced dots represent the centers of the space and time 
sample intervals for a straight wire, and the curves show the integra- 
tion paths for a given source segment on that wire as p is changed. 
Note that unless p = 0 and s coincides with a space sample point (i.e., 
the observation point is collocated with the sample point on the wire), 
the integration over a single space segment can involve more than one 
time sample of the current on that segment. This shows one reason 
for using basis functions that smoothly vary over the entire space-time 
domain. 


is limited to one or a few segments of length A, centered at 
observation point 7; and is equivalent to an applied voltage 
Viz = -4;8;;. The computation difference between a scattering 
and radiation calculation is minor. 

The effect of lumped or distributed loads is mathematically 
indistinguishable from the treatment of the exciting source. 
The load causes a voltage drop in opposition to (passive) or in 
the direction of (active) the exciting field. The actual load 
voltage depends on its V-J characteristics. For linear, passive 
loads we have (resistance R, inductance L, capacitance C) 


t 
Vioaa 7, =r, OR) +L(r) - U(r, t)+C(r) J q(r, 1) dr 


the effects of which can be readily included in the numerical 
treatment previously outlined. For example, this would result, 
in the case of our prototype integral equation, in 


8(x, t) > 8, 1)- Vioaa(x, t) 


and considering a resistive load only at x,, equation (6.5) 
becomes 


N 
By~ eZ, Ar ii | 
i'=1 , 
Ar Zin *RL (10) 
while all other A;;’, remain as given by (10). 

When the load is nonlinear, if for example the value of Rz, 
depends on J;,, then formally (10)’ still applies. But since 
R,, then depends on Az;, which in turn depends on R;, we 
generally must solve the equation by iteration at each time 
step. For the special case of an ideal diode, however, which is 
specified by only a forward and reverse resistance, we can 
readily solve (10)' by finding Az; for Ry = 0, then using the 
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direction of current thus determined to establish R;,, since 
the load cannot in this case reverse the current flow. We might 
also consider time-varying loads and other more general non- 
linearities [38]. 


B. Transform Techniques for Obtaining Transient Results 


Extensive literature on the numerical treatment of frequency- 
domain integral equations already exists [39], [40]. We, 
therefore, will not consider that topic any further here. 


C. Direct and Transform Techniques Compared 


In both the frequency domain and time domain, a system of 
equations is developed for sampled values of the unknown. 
The numerical solutions of these equations are substantially 
different, however, due to fundamental differences in their 
respective formulations. In the integral-equation treatment, 
for example, the interactions between the unknown sampled 
values are global in the frequency domain; i.e., the total field 
at a given observation point is due to the unknown sources 
distributed over the entire boundary. The spatial separation 
between the source and its field is manifested by a geometric 
attenuation and phase shift which are distance dependent. 
Consequently, all the unknown samples are mutually depen- 
dent, and must be solved simultaneously. The solutions are 
usually accomplished via matrix factorization or inversion. 

In a time-domain integral-equation treatment, on the other 
hand, the interactions between unknown samples are displaced 
in time by an amount equal to that required for a field to 
propagate between them at the speed of light. This displace- 
ment (time retardation) means that a particular unknown 
sample value at a given point in space and time is essentially 
determined by the exciting field at that same space-time point 
and by the scattered fields there which arise from earlier times 
and more distant locations. Consequently, the unknown sam- 
ples can be solved at any time step, provided all sample values 
at earlier times are already known. The time-domain problem 
is thus solved via time-stepping and without matrix inversion,* 
given the initial values of the unknown samples. 

Table I summarizes the main features of the direct time- 
domain and transform frequency-domain integral equations 
used to obtain transient results. The main differences between 
the two approaches are that: 1) the frequency domain provides 
source-independent results at a single frequency, so many fre- 
quency samples must be computed and transformed to find a 
transient result; 2) the time domain provides source-dependent 
results over the equivalent of a band of frequencies to give 
a transient result directly, but it must be repeated for each 
source; e.g., in computing the monostatic radar cross section. 
The corresponding computer times associated with their appli- 
cation are approximately [1] 


Ty (AzN?Ny + BINNTNg)Ny (11a) 


Ty © (AgN* + ByN? + CyNNy+ DgNNjNg)Ne (116) 


where A;,°"° , Dy are computer- and algorithm-dependent 
timing coefficients, and the subscripts t and f denote time- and 
frequency-domain quantities, respectively. The two terms in 
T, are due to current computation and far-field evaluation, 
where Ny is the number of time steps, NV, is the number of 


3 An initial matrix inversion may be necessary if interaction within 
the same time step of two or more space samples is allowed. 


1407 


TABLE I 
Frequency Domain Time Domain 
ane Maxwell's 
7 iw equations . Time dependent 
L(w)f(w) = g(a) Plus 8C, L(t) f(z) = g(t) 
etc. tT=t-R/e 
N N 
2 ace 
Zh = Gye ON 2. fay =O 
ik 7 
> k = 9 poply MoM 2 ikikj? = 944 
: to get 7 
i=1,...,N Nth order system Lel,...,N 
k= Ly..eyNy 
N N 
f, = Yy,Q.3 aN? Matrix foie 7 any 
- 2 1K? manipulation 4j 2 KG je? T 
yields - * 


Solution obtained for 
many sources but 
single frequency. 


Solution obtained for 
single source but for 
many frequencies. 


Do for 2 = 1,...,Ne 
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fis 


incident fields, and N, is the number of far-field evaluations. 
In Ty we similarly have terms in order of their appearance due 
to impedance matrix calculation, matrix factorization or inver- 
sion, current computation, and far-field evaluation, where Nr 
is the number of frequency steps. 

Upon introducing numerical values for the timing coefficients 
(for a CDC-6600), and relating the various numbers of current 
and field evaluations to the object size in wavelengths at the 
highest frequency of a transient calculation [1], we find the 
following computer-time approximations for wire (subscript w) 
and surface (subscript s) objects excited by a single source 
(Nz = 1): 


Tt,w ~ 10°'(LAmin)® s 

T3,5 ~~ 7(C/Amin)® § 

Tyw ~ (1+ 0.015(L/Amin)I(L/Amin)® 8 
Ty,5 ~ [2 + 0.83(C/Amin)? }(C/Amin)® 8. 


Here L is the total length of the wire object, C is the circum- 
ference of the smallest sphere which can contain the surface 
object, and Amin is the minimum wavelength for which the 
calculation is valid. Only the dominant terms from (11) have 
been retained. These results, plotted in Fig. 10, demonstrate 
the computational advantage of the time-domain approach to 
obtain transient results for a wire. For surface objects, on 
the other hand, the time-domain approach becomes more 
efficient only for C/Amin values 2 3. These results indicate 
general trends only, since they are based on nonoptimized 
codes and do not include the effects of symmetry or other 
influencing factors. They are useful, however, in drawing 
attention to an important aspect of transient calculations. 
Note that, for surface objects, computer time has a higher 
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Fig. 10. Computer time estimates to obtain transient results for a single 
excitation for both wires and surfaces as a function of length and 
circumference in wave-lengths at the highest frequency of the compu- 
tation, and using both a direct time-domain solution and a trans- 
formed frequency-domain solution. A significant advantage for the 
wire time-domain technique is shown. For surfaces, however, the dif- 
ference is not as great, perhaps because the time-domain code used 
for this comparison [13] re-computes the interaction elements at 
every time step. This data should be regarded as approximate and 
only a guideline. Possible reductions due to symmetry were not in- 
cluded. Note that the surface-object time estimates are larger than 
given in [1] because more conservative sample-density estimates were 
used here. 


order dependence on size because it involves an area, rather 
than a line, sampling of unknowns. 

A further computer-related characteristic that differentiates 
a time-domain from a frequency-domain approach for transient 
computation is the computer storage required. Because a given 
surface current and charge sample are needed during the com- 
putation only for the time when their fields interact with the 
test of the object, the storage required can be considerably 
less than the NN 7 otherwise required [1]. 


D. Modeling Guidelines 


Because time-domain techniques involve more user decisions 
than their frequency-domain counterparts, in terms of tem- 
poral sampling and the time variation of the exciting field, 
they can be somewhat more demanding to use. If, however, 
the goal of the frequency-domain approach is to develop tran- 
sient behavior via a transform, then similar decisions regarding 
frequency sampling and the spectrum of the exciting field are 
required. In each approach, these factors are important to 
both the efficiency and accuracy of the computation, but it 
temains unclear how to choose them a priori to most nearly 
attain the conflicting goals of computational speed and numer- 
ical accuracy. Also common to each is the need to spatially 
sample the object being modeled in a manner consistent with 
the temporal sampling and the accuracy-efficiency aspects of 
the calculation. 

In our discussion of sampling questions above, the sampling 
density used the Shannon-Kotelnikov (also sometimes called 
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Nyquist) rate of two samples per period of the highest fre- 
quency component in the waveform (or spectrum). But this 
value is an upper bound on the temporal-sampling interval and 
is attainable only when the basis function exactly matches the 
components in the waveform. In our example, the exponential 
time variation discussed in connection with SEM results in a 
damped, time-harmonic oscillation which exactly represents 
the late-time free response of a conducting object so that the 
Nyquist rate does apply. Thus the temporal- and frequency- 
sampling intervals have their maximum theoretically possible 
values when using the SEM representation. SEM also provides 
the most efficient representation and transform between the 
time and frequency domains for complex exponential func- 
tions via the Laplace transform pair. 

However, when the frequency and time domains are con- 
nected via the Fourier transform pair, the sampling intervals in 
time and frequency must be decreased because the basis func- 
tions (in this case, exponentials with purely imaginary argu- 
ments) are not the most appropriate for a waveform (or spec- 
trum) of complex exponentials. One can better understand 
this fact by observing that a complex exponential can only be 
approximated, but not exactly matched, by sinusoids unless 
an infinite spectrum of such sinusoids is used. In practice, 
when using the FFT or alternate numerical transforms, we find 
typically that 10-12 samples per period of the highest fre- 
quency component are required to obtain acceptable accuracy. 

The above discussion was included to emphasize that sam- 
pling densities are determined by several factors, and that 
theoretical expectations based on optimum conditions may 
not be met. Here we conclude our discussion by summarizing 
the following modeling guidelines concerning time-domain cal- 
culations for wires using the EFIF: 

1) space sample A should satisfy 


2a <AS Xmin/6 


where a is the wire radius and Amin the shortest wave- 
length of concern; 
2) time sample 6 should satisfy 


6 S Amin/6c 


where a lower limit is not given because the theoretical 
restriction that applies to the magnetic-field case where 
& 2 A/c does not apply here [15]. 

These sample intervals are smaller than the Nyquist rate, be- 
cause the space-time basis function used (a nine-term second- 
order polynomial) is only an approximation of the actual 
current and charge variations being modeled. 


£. Degree of Applicability 


Problem types for which a computer model may prove suit- 
able depend on the user’s judgment of computer costs involved, 
the degree to which the model resembles the problem of 
interest, etc. The prospective user should be aware of various 
modeling aspects which affect potential applications. We 
categorize them as configurational factors (relating to the 
object), environmental factors (relating to the object’s sur- 
roundings), and miscellaneous factors. 

1) Configurational Factors: As a general rule, wire objects 
are more amenable to computer modeling than are surface 
objects, because they involve a linear rather than an area sam- 
pling. Wire objects successfully modeled in the time domain 
include straight and curved (modeled as piecewise linear) wires, 
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open- and closed-ended wires, wires having bends and multiple 
junctions, and wire grids. Surface objects modeled include 
smooth closed surfaces, closed surfaces with edges and vertices, 
and open surfaces (shells). Hybrid objects (i.e., objects having 
wire and surface features) have also been modeled. 

2) Environmental Factors: The characteristics of the me- 
dium in which an object is located can have a very great impact 
on the feasibility of developing a computer model. Most 
results obtained to date apply to infinite, uniform, lossless, 
nondispersive, linear, and isotropic media. A relatively simple 
extension is possible to handle up to three orthogonal, per- 
fectly conducting image planes. Some work has also been 
done for a general, nonlinear medium [41]. 

3) Miscellaneous Factors: Other time-domain features avail- 
able include the ability to handle lumped and distributed 
resistance, inductance and capacitance, and nonlinear loads. 


F. Limitations 


The numerical implementation of direct time-domain tech- 
niques can be generally inferred from the list of applications 
mentioned above. But listing a given problem does not neces- 
sarily mean that it can be modeled without difficulty. Many 
problems described are one of a kind and include features for 
which the treatment remains uncertain (e.g., sharp bends on 
wires, edges and vertices on surfaces). Furthermore, the solu- 
tions which have been worked out may be for special cases and 
not yet generalized (e.g., expansion for an axi-symmetric 
surface object for axial, but not oblique, incidence). 

Perhaps the most serious limitation concerning direct time- 
domain computations, however, is the computing capability 
required. Even the CDC 7600* computer can be challenged 
by apparently simple wire calculations. For example, with a 
core of ~5 X 10° words available for variable storage, the time 
response of a straight wire having more than 150 space samples 
cannot be carried out in core for the minimum of four transit 
times required to obtain its poles (assuming A=c6). Looked 
at another way,> 2NN7 + LON? < 5 X 105, so that with Nr ~ 
4L/c5 =4L/A and N=L/A, we have A26X10°L. If 
samples are saved for only the time span L/c for which they 
are needed, then A can be reduced by only a small factor of 
less than 2, allowing a slight increase in the equivalent fre- 
quency range. If only the peak response or the energy col- 
lected is sought, then the computation might be completed 
much earlier than L/c. Generally, however, wide-band calcu- 
lations can require substantial storage (and time) resources. 


VII. EXPERIMENTAL ASPECTS OF DIRECT TECHNIQUES 


A presentation of transient electromagnetics would be in- 
complete without mentioning the role of experimental mea- 
surements. Of course, as is true for analysis and computations, 
we can derive transient results from either direct or transformed 
measured data. Our concern here is with the former since 
frequency-domain experimentation is well known. In the 
discussion that follows, we outline a general approach to time- 
domain measurements, describe one particular facility, present 
some representative results, and summarize various applications 


“Reference to a company or product name does not imply approval 
or recommendation of the product by the University of California or 
the U.S. Energy Research and Development Administration to the 
exclusion of others that may be suitable. 

‘The 2NN7 arises from storing current and charge information, and 
the 10N? is due to other storage requirements. 
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for which such techniques might be considered. This approach 
was pioneered at Sperry Research Center [33], [42]. 


A. Time-Domain Measurements 


The most significant difference between time- and frequency- 
domain measurements is that the former uses time-range gating 
to eliminate the need for the reflection-free environment the 
latter requires, A time-domain measurement is conceptually 
quite simple, requiring as a minimum an impulsive electrical 
source, a sensor, wide-band sampling oscilloscope to measure 
the sensor output, and associated cabling. Source characteris- 
tics primarily determine the overall bandwidth and signal-to- 
noise ratio of the transient experiment. Commercially avail- 
able pulsers having rise times approximately 250-300 ps and 
peak voltages approximately 1-2 kV give useful maximum fre- 
quencies extending to 3-4 GHz. Pulsers with even faster rise 
times, approximately 100 ps, are available, but they have lower 
outputs and consequently lower signal-to-noise ratios. 

The sensor used for time-domain measurements can be simple 
or elaborate, depending on the application. For determining 
the transfer or self-admittance of a wire antenna, for example, 
the antenna itself serves as the sensor. For determining surface 
current or charge density, we can use various probes based on 
loops or monopoles. If these probes are small relative to the 
object tested, their outputs will approximate the time deriva- 
tive of the local magnetic or electric field. Therefore, they 
are often called B-dot and E-dot sensors, respectively. The 
Air Force Weapons Laboratory recently developed a line of 
precision-designed and constructed sensors in lead arrange- 
ments for various applications. 

Available sampling oscilloscopes can measure to 12-14 GHz 
and thus cover the bandwidths provided by current and antici- 
pated pulsers. Single pulse transient digitizers are also avail- 
able, but their equivalent bandwidths are generally less (e.g., 
1 to 2 GHz), the number of data samples they can provide is 
more limited and their costs are greater. 

Although only the items described above are essential for 
time-domain measurement, other components can crucially 
improve the accuracy and/or efficiency of a system. Perhaps 
most important is a minicomputer for data acquisition and 
processing. It can monitor and control the experiment and 
perform initial data processing such as smoothing, averaging, 
or FFT-transforming to the frequency domain. Coupled with 
peripheral equipment such as a magnetic tape unit or paper 
tape reader/puncher, it permits data storage for post processing. 
The data acquisition system can be supplemented by an X-Y 
plotter for real-time data presentation in hard copy. 


B. The LLNL Transient Facility 


Fig. 11 is a block diagram of the Lawrence Livermore Na- 
tional Laboratory (LLNL) transient system [34], including all 
the electronics and a schematic depiction of the ground-plane 
range, transmitting antenna, and test target in a typical mea- 
surement configuration. An advantage of the ground-plane ar- 
rangement is that the instrumentation cables can be introduced 
from under the ground plane, and are thus kept from interfer- 
ing with the measurement. A disadvantage is that only objects 
having a symmetry plane can be measured, and polarization 
and incidence angles are limited. Time-domain measurements 
can be made in a free-space configuration as well, but, because 
isolating the cables is difficult, fiber optics or other nonmetallic 
data links may be necessary. Otherwise, only the early-time 
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Fig. 11. Block diagram of a transient, ground-plane range shows its 
conceptual simplicity [33], [34]. Major elements are the ground 
plane, pulser, signal-sampling oscilloscope, and computer data- 
acquisition, processing, and control components. 


portion of the measured waveform may be uncontaminated by 
cable-induced artifacts. 

The range clear time, which determines the low-frequency 
cutoff of the measurement, is established by the first arrival at 
the measurement point of a reflection from the walls or ceiling. 
In the LLNL range, the clear time is 20-25 ns, although for 
some measurements useful data can be obtained over time inter- 
vals of 200 ns. Several other factors affect the overall quality 
of a given measurement, including the degree of wavefront 
planarity achievable over the target, the target-size-to-wave- 
length ratio at the highest effective frequency, the signal-to- 
noise ratio established by the pulser output and instrumenta- 
tion, pulser time and amplitude stability, sensor sensitivity, and 
target-response characteristics. These factors are discussed 
more completely elsewhere [34]. 


C. Some Representative Results 


Results obtained from two kinds of measurements performed 
on the LLNL range are compared in Figs. 12 and 13 with cal- 
culations based on a thin-wire time-domain computer code. 
Fig. 12 shows the input resistance and reactance of a base-fed 
monopole as functions of frequency obtained via an FFT from 
time-domain data. The experiment was performed by measur- 
ing the pulse reflected from a cable terminated at the ground 
plane with and without the monopole plugged into it, and 
computing the frequency-domain reflection coefficient from 
which the impedance is obtained. 

The data in Fig. 13 are the transient response of a V-loaded 
monopole, presented as the current through a 50-22. load. The 
numerical model includes the transmitting antenna (in this 
case a long (~4m) straight wire), for which the numerical 
exciting voltage is obtained from measuring the input volt- 
age to it. Apart from a slight difference in the time base be- 
tween the computed and measured result, they are in excellent 
agreement. 
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Fig. 12. Monopole impedance [(a) resistance and (b) reactance] as a 
function of frequency as obtained from a time-domain calculation 
(—) and from transient-range measurement (***) transformed to the 
frequency domain using a FFT (29]. 
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Fig. 13. Current induced on a V-loaded dipole as computed (***) using 
a time-domain approach and as measured (——) on the transient range 
{34]. The calculation model includes the transmit antenna (a long 
monopole) and the measured voltage used to excite it. No normaliza- 
tion of the two results was performed; their absolute values are shown. 


Use of a time-domain measurement to estimate the error 
associated with thin-wire modeling is demonstrated in Figs. 14 
and 15. There we show three models for a 747 aircraft. One 
is a scale model, and the other two are circular-cylinder models, 
one having a diameter comparable to the scale model (the pipe 
model) and the other having a smaller diameter (the thin-wire 
model). Each model was mounted nose-down, perpendicular 
to the ground plane, and the current excited at its connection 
point by an impulsive wave was measured with the results 
shown in Fig. 15(a). The corresponding spectral responses 
obtained by Fourier transforming the time-domain data are 
presented in Fig. 15(b). 
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Fig. 14. Physical and numerical modeling errors can be assessed using 


models like these. 


Shown here are a scale-model and pipe and thin- 


wire models of a Boeing 747 aircraft [34]. 
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Fig. 1S. Response of models in Fig. 14 can be used to evaluate the 
physical modeling error. 


These results are useful to evaluate both the physical model- 
ing error €, which arise from replacing the physical problem 
with an idealized mathematical representation, and the numer- 
ical modeling error ey which arises from obtaining an approxi- 
mate solution for the numerical model selected. Comparison 
of the pipe-model and thin-wire model results with those for 
the scale model provide information concerning €, (Fig. 15), 
and comparison of thin-wire measured and computed data 
(Fig. 16) provides information concerning €y. While equivalent 
data could be obtained from frequency-domain measurements 
and computations, those from the time-domain are generally 
more appropriate and efficient for transient applications. 


—— Experimental 
--= Computed 


Nose current ~— arbitrary units 


0 4 8 12 16 20 
Time — ns 
Fig. 16. The numerical modeling error can be examined by comparing 
measured and computed results for the numerical model. In general, 


as observed here for the thin-wire model of the 747, the numerical 
error can be made less than the physical error. 


D. Applications 


Two kinds of direct time-domain measurements (input im- 
pedance and transfer admittance) are illustrated above. The 
range of potential or actual measurement types is quite broad, 
and can be divided into two main categories, closed and open 
systems. 

I} Closed Systems: Following the lines of time-domain 
teflectometry, we can measure constitutive parameters, system 
tesponse functions, etc. [33]. In closed geometries, such as 
coaxial lines and waveguides, geometric attenuation of the 
exciting field is not the problem it can be in open systems. 
Therefore, we can attain wider measurement bandwidths in 
the closed system with a given pulse spectrum. 

2} Open Systems: The range of applications available in an 
open system includes transfer function, scattering, and input- 
impedance measurements. Input-impedance measurements 
are probably least demanding of these, because the input and 
output signals are not geometrically attenuated. The transfer- 
function measurement can be used for obtaining the trans- 
mitting and receiving patterns of antennas, while scattering 
Measurements can be monostatic or bistatic. Measurements 
where the exciting field is incident on the sensor may reduce 
the dynamic range available to observe the scattered field, 
because the sampling oscilloscope must be set on a range deter- 
mined by the maximum input signal. 
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The applications of time-domain measurements are of course 
similar in many ways to their frequency-domain counterparts. 
They are useful to validate analytical results and computer 
codes in either the time domain or frequency domain. More 
important, they can be used to generate data not derivable 
from analysis and thus provide information not otherwise 
available. 


VIII. REPRESENTATIVE RESULTS 


In the preceding discussion, we have briefly summarized 
various physical mathematical, numerical and experimental 
aspects of direct, time-domain modeling. We continue our 
presentation by including below representative numerical 
examples to demonstrate features essentially unique to time- 
domain treatments, and to show various types of applications. 


A, Features Unique to Direct Time-Domain Analysis 


1) Nonlinear Capability: Solution directly in the time do- 
main permits analysis of nonlinear problems by a time-stepping 
procedure. This approach has been used with finite-difference 
techniques to handle the problem of a metallic object im- 
mersed in a nonlinear medium [41]. Schuman [42] and others 
[44], [45] have considered a thin straight wire loaded at its 
center by a diode. This case, as well as nonlinear loading of 
more general wire geometries, is possible using the approach 
developed by Miller et al. [17], which is the one we concen- 
trate on here. 

Figs. 17 and 18 show the response of a linear dipole antenna 
loaded with a diode. The diode was placed in series with a 
Gaussian voltage source. Fig. 17 shows early-time history of 
the current and charge distribution along the structure. The 
diode is initially biased in the forward-conducting direction, 
allowing a charge separation to build up across it. When the 
exciting voltage decreases to zero, the accumulated charge 
then reverse-biases the diode, allowing only a small leakage 
current in the reverse direction. The dipole then responds 
as two shorter dipoles placed end-to-end but insulated from 
each other. Comparison of the center current for this case 
(Fig. 18(a)) and the unloaded case (Fig. 2) illustrates that the 
initial current is essentially the same for both (there is a slight 
amplitude difference caused by different wire radii) until the 
current in the linear case changes sign. From then on, the 
two currents are distinctly different. The late-time far field 
tadiated broadside from the dipole is predominantly a damp- 
ened sinusoid at twice the frequency of the unloaded dipole, 
as shown in Fig. 18(b). 

Two other examples of the capabilities provided by time- 
domain analysis for modeling nonlinear problems are demon- 
strated in Figs. 19 and 20. In Fig. 19 are presented the input 
current and broadside radiated field for a dipole antenna 
which is loaded with diodes, one for each of the segments used 
to model it. The current naturally flows in only one direction 
as can be seen in Fig. 19(a), and the resulting radiated field, 
Figs. 19(b) and (c), is an extended doublet. The first pulse is 
due to the initial turning on of the current, and the second is 
due to stopping of the current as it reaches the end of the 
antenna. Because of the diodes, the flow of current (and 
charge) in the opposite direction is stopped so that no further 
pulses of radiation are produced. By adjusting the diode 
characteristics, their number and locations, the antenna length, 
the exciting voltage, etc., some control of the radiated field is 
possible. 

The result of illuminating by a plane wave a straight wire 
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Fig. 17. Nonlinear loads can be easily handled in the time domain. Here 
the current (—) and charge (—---—) distributions on a 1-m center-fed 
dipole having a diode load in series with the generator are shown at 
several instants of time. Observe that the diode presents zero imped- 
ance only to current flowing in the initial direction, and that it effec- 
tively acts as an open circuit in the opposite direction (impedance = 
100 k&). The two halves of the dipole thus retain a static charge dis- 
tribution at late times. 


> 


2 


Source current — mA 


Radiated fields times distance — V 


Fig. 18. (a) The transient feedpoint current on the nonlinearly loaded 
antenna of Fig. 17. (b) The radiated field possesses a fundamental 
frequency twice that of an unloaded dipole. 


loaded at its center by a resistor whose resistance varies sinu- 
soidally in time is demonstrated in Fig. 20. Shown there are 
the time-dependent load current from the time the incident 
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Fig. 19. Potential for pulse shaping, using nonlinear distributed loading. 
A 1-m dipole, modeled numerically by 60 segments, is loaded with a 
series of ideal diodes, one per segment. (a) Application of a Gaussian 
voltage pulse across the center segment produces this feedpoint cur- 
rent. (b) Broadside radiated field. (c) Frequency spectrum. Note the 
sharply curtailed time response of the field caused by the diodes and 
the absence of late-time ringing. 


wave is turned on and the scattered field, as well as the Fourier 
transform of the latter. The effect of the load modulation 
caused by the resistor is to produce scattered energy in side- 
bands of the incident-wave frequency. Depending on the load 
properties available, it may be possible to transform the inci- 
dent field frequency to selected sidebands, which meet speci- 
fied characteristics, e.g., to simulate a moving target via a simu- 
lated Doppler-shifted radar return. 

2) Time Gating: Time gating can be used in transient EM 
measurements to eliminate the effects of reflections from 
scatterers physically separated from the target. This powerful 
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Fig. 20. Time varying load in a time-domain computation, illustrated 
by a 9-m dipole that is illuminated by a broadside incidence of a 
16-MHz-modulated Gaussian pulse. Load resistance at the center is 
square-wave modulated. Scattered energy appears in bands, harmoni- 
cally related to the modulation and illumination frequency. (a) Load 
current versus time. (b) Scattered field versus time. (c) Scattered 
field magnitude versus frequency. This shows the possibility of scat- 
tering signals in desired bands other than the frequency band of the 
incident signal. 


technique makes possible wide-band measurements without 
expensive anechoic chambers, etc. But time gating can also 
be used numerically, as demonstrated in Fig. 21 [46]. Here, 
a dipole of finite length is illuminated by the time derivative 
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Fig. 21. Time-range gating can be used to separate reflection effects 
from other phenomena. (a) The computed current prior to end reflec- 
tion at the center of a long (120-m) center-fed dipole modeled with 
120 segments is extrapolated in time to approximate the current that 
would be seen on an infinite antenna. (b) The transfer function can 
be obtained by transforming the current to the frequency domain. 
The transfer function compares well with analytical results. Several 
calculations were required to span the ka range shown [46]. 


of a Gaussian plane wave at broadside incidence. Fig. 21(a) 
shows the current at the center of the dipole. The calculation 
ended before end reflections arrived. Thus time gating permits 
us to obtain an infinite-wire response from a finite-wire re- 
sponse. Transformation to the frequency domain establishes 
the validity of this approach for the transfer admittance. Fig. 
21(b) shows the analytic results and the FFT of Fig. 21(a). 
We can find the frequency range over which this technique 
produces valid data in the same way as in transient measure- 
ments. The clear time determines the lowest frequency of 
valid data, and the sampling intervals determine the highest 
frequency. 

Besides providing a way to circumvent the effects of un- 
wanted reflections, time gating can also be used for detailed 
study of such reflections. Consider as one example, the 
center-excited wire antenna whose space-time current be- 
havior is presented in Fig. 1. The spatial variation of the cur- 
rent and charge pulses which are incident on the wire’s ends 
can be separated in time from their reflected counterparts, 
to obtain the frequency dependence of the end-reflection 
coefficient. 

The equivalent information could also be obtained from 
the time variation of these same pulses at a given location on 
the wire, an example of which is shown for the center-excited 
straight-wire antenna in Fig. 2. A more interesting case is illus- 
trated in Fig. 22(a), which displays the time variation of the 
current at the center feedpoint of a bent-wire antenna, formed 
by bending a straight wire at right angles one quarter of its 
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Fig. 22. A useful feature of time-domain data is the possibility it offers 
for separating in time the effects of various structure features on cur- 
rent and charge pulses. An example is shown in (a) for the current at 
the center of a wire which has right-angle bends one-quarter of its 
length from each end. In part (b), corresponding results are presented 
for a straight wire having 600-2 loads one quarter of its length from 
each end. Both the incident and reflected pulses are clearly demar- 
cated, so that their individual Fourier transforms could be obtained 
and used to develop the reflection coefficients for the bends and loads 
as a function of frequency. Because the bends electrically resemble 
a short and the resistors an open-circuit termination, the current 
pulses they reflect are of opposite sign. 


length from each end. Clearly observable are the separate 
pulses reflected from the bend and the end of the antenna. 
The effect of adding 600-Q resistors to a straight wire one- 
quarter of its length from the ends and excited at its center 
is shown in Fig. 22(b). Note that the reflected current pulse 
is of opposite sign from that of Fig. 22(a), because the impe- 
dance discontinuities in these two cases are similar to short- 
and open-circuit conditions, respectively. 

3) Physical Insight: Transient methods can be advanta- 
geously used to help understand the physical characteristics 
which underlie a structure’s EM response. The insight thus made 
possible may be greatly enhanced by novel data-presentation 
techniques. For example, a motion picture of the currents 
along the structure and the radiated fields has been found 
useful in.showing differences between radiation associated 
with end reflection compared with that from a curved wire. 
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Fig. 23. Space-time contour plots of the current can convey a great 
deal of information concerning transient behavior. The time variation 
is shown on the vertical axis and the space variation is shown on the 
horizontal axis, where the scales are in the ratio of c, so a 45° straight 
line represents motion at the speed of light. Results for the center- 
excited (a) straight-wire, (b) bent-wire, and (c) resistively loaded 
straight-wire antennas of Figs. 2 and 22 are shown here. The differ- 
ence in radiation mechanisms for these antennas, and other features 
as well, can be inferred from these results, as discussed in the text. 
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A motion picture format is not always suitable, of course, 
but the presentation continuity it provides in both space and 
time can be obtained in other ways. One example is presented 
in Fig. 23. There we present constant-current contours as a 
function of space and time for three thin-wire antennas. These 
include the center-excited straight wire of Fig. 1, the bent-wire 
antenna of Fig. 22(a), and the resistively loaded wire of Fig. 
22(b). Such results represent essentially complete current solu- 
tions from which current plots like those of Figs. 1 and 2(a) 
could’ be derived. On these plots, distance (horizontal axis) 
and time (vertical axis) are scaled such that a pulse moving at 
the velocity of light describes a straight line with a slope of 
unity. 

The straight-wire result (Fig. 23(a)) illustrates that the inci- 
dent current and charge pulses undergo little decay traveling 
along the wire, as shown by the closing of all the contours near 
the wire’s end, but that the reflected pulse is slightly less than 
that incident; energy has been lost to radiation. For the bent 
wire (Fig. 23(b)), the effect of the bends shows up in two ways, 
in a peaking of the current at the bend (note the closed con- 
tours) and in the reflected current pulse (note the contours 
sloping inward from the bend’s location). ‘The results for the 
resistively loaded straight wire (Fig. 23(c)) also exhibit reflec- 
tion from the load, accompanied by a decreased current ampli- 
tude there. In all cases, these current contour plots provide an 
alternative, more detailed view than can be obtained from plot- 
ting the current at a single point as presented in Figs. 2(a), 
22(a) and 22(b). 

Space-time contour plots depicting the far field are shown 
in Fig. 24 for the three thin-wire antennas of Fig. 23. The 
purpose of presenting the far fields in this manner is to provide 
completely their variation with angle and time in a fashion 
analogous to the current contour plots. Also, and of perhaps 
greater interest, is the opportunity such a format presents to 
demonstrate clearly where the radiation originates from the 
antenna. Note that with this format, an impulsive source 
located at the coordinate system origin will produce a circular 
ring showing it to be located at the center of a spherical shell 
of radiation. Any movement of a source away from the point 
from which the time is measured will produce a elliptical pat- 
tern of the kind indicated in the figure. 

The far-field contour plot for the straight, center-excited 
antenna (Fig. 24(a)) shows clearly the initial burst of radiation 
from the source region, and the subsequent radiation from the 
ends. For the bent-wire antenna, additional far-field contribu- 
tions result (Fig. 24(b)), showing that radiation occurs due to 
teflection from the bends. A similar result is found for the 
resistively loaded straight-wire antenna (Fig. 24(c)). These 
plots identify explicitly the sources of far-field radiation and 
show the interesting possibilities made available from time- 
domain solutions. 

4) Broad-Band Nature of Time-Domain Results: The fre- 
quecy-domain analysis of a structure is naturally limited to 
discrete-frequency samples. The user of a code must choose 
a priori those frequencies where the structure is anticipated 
to have interesting or useful behavior. An incorrect choice 
results in the need to continue the calculations until the desired 
information is obtained, and confidence has been acquired 
that essential aspects of the object’s behavior have not been 
overlooked. This is especially the case when the narrow-band 
resonances of a high-Q structure are being studied. 

A single time-domain analysis, on the other hand, is inher- 
ently broad band. As long as the excitation contains sufficient 
energy in the frequency range of interest, the entire response 
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Fig. 24. Examination of the far field can reveal the source locations on an object. Here we present contours of constant far-field values as a 
function of angle and time (or space) relative to the feedpoint of a wire antenna excited at its center by an impulsive source. In (a) is shown 
the result for the unloaded straight wire, where we deduce that the initial, circular band is due to radiation from the source, with the inner, 
nearly circular bands caused by radiation from the ends. The radiation contribution due to right angle bends located at the middle of each 
half of the dipole is demonstrated in (b). In (c) is shown the effect of adding resistive loads midway between the antenna’s center and ends. 
Both these results are similar to those of (a) except that a new radiation band occurs at t ~ L/4c, obviously caused by the bends and loads, 
respectively. For all cases shown, the antenna lies on the horizontal axis of the figure with its center at the vertical mark, and it is the @ com- 
ponent of the field which is plotted. In (b), the bent ends of the antenna are normal to the paper. Note 0—observation angle. r—observation 


distance. t—observation time. 


over that range is automatically obtained. For example, con- 
sider a frequency-independent antenna, the conical spiral. This 
frequency independence occurs, however, over some limited 
bandwidth where the antenna was designed to operate. If one 
needs characteristics over a wider band (for pulse applications 
or EMP coupling analysis), transient analysis conveniently 
provides them. In Fig. 25 we show some characteristics of a 
typical conical spiral. Fig. 25(a) depicts the feed-point current 
resulting from a Gaussian excitation in series with a 350-Q 
source impedance [47], included to damp out the low-frequency 
resonances that would otherwise cause the response to persist 
for many transit times. A current contour plot for this antenna 


is shown in Fig. 25(b). Note the closed current contours due 
to radiation damping and increasing pulse width caused by 
dispersion, as the current pulse propagates along the curved 
wire of this antenna. Transforming the current of part (a) to 
the frequency domain yields the input impedance (corrected 
to remove the generator resistance) presented in Fig. 25(c). 
The calculation produces low-frequency resonances of very 
high Q. The frequency-independent region of this particular. 
conical spiral is fairly narrow, from 1.2 GHz to about 3 GHz. 
At low frequencies, the spiral acts very much like an open- 
circuit transmission line because it is basically a twisted pair 
of wires, with the wire separation increasing slowly. Conse- 
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Fig. 25. (a) Feedpoint current for a conical-spiral antenna excited by a Gaussian pulse. The effect of a 350-2 series resistor at the feed point 
is included. The numerical model consists of a piecewise linear approximation (60 segments) to the actual antenna, whose total length is 1 m 
and radius is 1 mm, and which has two arms of 2.5 turns wound on a cone of half-angle 10° with small-end and large-end diameters of 3.36 cm 
and 9.16 cm, respectively. (b) Space-time contour plot of the current on this conical antenna, demonstrating radiation damping of the cur- 
rent (closing of the current contours) and dispersion (broadening of the pulse). (c) Impedance for a conical-spiral antenna excited by a Gauss- 
ian pulse, obtained by transformation of (a) to the frequency domain. The late-time ringing of the spiral that would otherwise occur is 
damped out by the resistor and has been subtracted from the resistance curve. 


quently, the observed resonances are expected and qualitatively 
explainable from transmission-line theory. 

5) Use of Time-Domain Data for Object-Pole Finding: A 
convenient procedure for finding the EM poles of an object is 
to apply Prony’s method to some measure (current, charge, 
near field, far field) of its transient response. One especially 
useful characteristic of these poles is their independence of 
the excitation, a subject demonstrated elsewhere [29], [30]. 
Here, we demonstrate some other interesting aspects of EM 
poles. 

Figs. 26 and 27 show results of resistively loading a dipole 
at its center, and bending it at a right angle at its center. The 
transient waveforms are the backscattered fields from the di- 
poles when illuminated by a Gaussian-pulse field. In the former 
case, the a = 1, 3, 5,’ °° poles are more lossy than the a = 
2, 4, 6,° °° poles, while the converse is true in the latter 
case. This occurs because the odd-numbered modes are even 
about the dipole’s center, so that the resulting current maxima 
at the resistive loads produce a large dissipative loss, while the 
even-numbered modes are odd and have charge maxima at the 
center, thus producing a large radiative loss at the bend. 

Besides being useful to find the pole locations for SEM, 


Prony’s method serves as well to store the transient waveform 
in a shorthand form (Fig. 28). One interesting application of 
data handling is the ability to extrapolate both forward and 
backward in time. Fig. 28(a) shows a portion of the computed 
backscattered field observed from the wire object shown. Fig. 
28(b) plots the pole locations obtained using only that portion 
of the response indicated in Fig. 28(a). The poles from Fig. 
28(b) (and residues which are not shown) were then used to 
fill in the entire transient response. This extrapolated response 
is compared to the original computed response in Fig. 28(c), 
with excellent agreement. 

6) Diagnostic Property of Time-Divergent Solutions: Many 
pitfalls may be encountered in applying any numerical tech- 
nique. A peculiarity of direct time analysis that can be used to 
advantage is the way the results display ill-conditioning. For 
example, experience shows that using the thin wire approxi- 
mation (see Miller et al. [17]) that wire segments should be at 
least as long as their diameter. In Fig. 29 we show the current 
on a dipole excited by a broadside-incident, Gaussian-pulse 
plane wave for two wire radii. In each case, the dipole length 
is 1 m, and the segment length (A) is 1.666 cm. The two curves 
in Fig. 29 show the result of increasing the wire radius, by 
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Fig. 26. Energy loss mechanisms can affect the pole locations in the 
complex frequency plane.”” Here we present the pole sets for a 60-m 
wire, showing the result when the wire is straight and loaded with a 
100-Q resistance at its center. 


only ~5 percent with A2a ~ unity. For the larger A/2a, the 
current is well behaved and not a function of the number of 
segments. (Of course, if too few segments are used, the results 
will be inaccurate.) But an exponentially growing current is 
obtained below some critical value for A < 2a, signifying that 
something has become ill-conditioned. The particular source 
of this phenomenon in the numerical procedure has not been 
found, but is apparently a result of the thin-wire approxima- 
tion. While the impedance matrix or interaction coefficients 
display nothing obviously wrong, the growing currents clearly 
signal that something in the numerical model is in error, how- 
ever. Similar instabilities have been observed if only a single 
segment is too short, or if the match point on one segment 
lies inside the volume of another segment, and thus serve as a 
diagnostic of an invalid numerical model. 


B. Applications 

1) Transient-Pulse Shaping: For some applications, it may 
be desired to radiate a specific transient-pulse shape. Physical 
constraints may prohibit the antenna from being frequency 
independent throughout the band of the desired pulse, and the 
antenna then functions as part of the wave-shaping network. 
For example, consider a zero-impedance voltage source con- 
nected to a conical spiral antenna (the same antenna used for 
the results in Fig. 25) [47]. Assume that we want the radiated 
field in one polarization plane to be the third derivative of the 
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Fig. 27. Here, the wire is unloaded but has a 90° bend at its center. 
The a=1,3,--+ poles are more lossy for the loaded wire, due to 


dissipative loss of those modes which have current maxima at its cen- 
ter, and the a=2,4,-+-- modes are more lossy for the bent wire, 
because the charge maxima at the bend produce a greater radiation 
loss. 


Gaussian, i.e., 
a 
Eraa(t)@ ie [exp (-a71?)]. 


The desired quantity is the voltage required to be applied to 
the antenna to produce that far-field behavior. To find this 
quantity, we first calculate the transfer function between the 
radiated field and the applied voltage, which is shown in Fig. 
30(a) for the conical antenna previously considered. Next, 
we obtain the spectrum of the required voltage by dividing the 
spectrum of the desired radiated field by this transfer function. 
The spectrum of the desired field must go to zero with de- 
creasing frequency faster than the transfer function, because 
the antenna cannot radiate a static field. We then find the 
required transient voltage (Fig. 30(b)) by transforming this 
spectrum to the time domain. To verify the result, we use this 
waveform in the time-domain wire model to obtain the radi- 
ated field shown in Fig. 30(c) from the direct time-domain 
calculation. This waveform is the desired third-time derivative 
of the Gaussian as specified. A similar approach could also 
be used with measured antenna characteristics. This kind of 
synthesis could also be performed using the frequency-domain 
of course; equivalent information is needed in either case. 

2) Wire-Grid Models: A reasonably complicated grid prob- 
lem is depicted in Fig. 31(a) of a wire-grid model of a light 
truck with a rear-mounted 108-in whip antenna. The effect of 
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Fig. 28. Extrapolation of a transient waveform can also be accom- 
plished using the poles extracted from it. (a) Part of a computed field 
scattered from the thin-wire object is illustrated. (b) Poles extracted 
from a portion of the transient waveform. (c) The waveform recon- 
struction, which compares well with the original waveform both be- 
fore and after the sampled portion. 


the ground was included as a perfect image plane. Fig. 31(b) 
shows the results for the input admittance obtained from a 
time-domain calculation. An input resistance of ~70 Q occurs 
at 27 MHz, the operating frequency of the transmitter, which 
may be recognized as the Citizen’s Band (CB) frequency. This 
result has been verified by actual operation, demonstrating 
the utility of the wire-grid model. 

A variety of other wire-grid calcuiations have also been per- 
formed using the time-domain approach. These include mesh 
models for plates, a fan-type antenna, and conical shell. Be- 
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Fig. 29. Solution instability caused by use of too-short segments. 
(a) The transient current on a 1 m straight wire modeled using 60 seg- 
ments (A = 1.66 cm) is displayed here for a wire 9 mm in radius, so 
that A/2a= 0.926. (b) Transient current for a wire 9.5 mm in radius 
where A/2a = 0.877. The divergent solution that occurs is evidently 
due to numerical inaccuracy of the thin-wire approximation that 
occurs with decreasing segment-length-to-radius ratio. 


cause wire grids are only physical approximations to closed 
surfaces, care must be exercised in their use, however.. 

3) Transient Response: Another obvious application of 
time-domain analysis is the direct (versus using a transform 
technique) calculation of coupling such as that due to EMP 
which of course is a transient phenomenon. For example, in 
Fig. 32 we show the response of a 10-m straight wire resistively 
loaded at its center illuminated from broadside by an EMP 
pulse. 

The ease of obtaining these responses conveniently permits 
efficient parameter studies. Fig. 33 shows the results of one 
such study. Here, wires of different lengths were exposed 
to an EMP pulse. The wires were loaded resistively at their 
centers, and the value of resistance was also considered a pa- 
rameter. For each combination of length and resistance, the 
transient calculation was performed and the peak load current 
and total energy dissipated in the load were evaluated. The 
results of these calculations are shown as contour plots, where 
the loci describe parameter pairs that produce the same peak 
current or total energy. Such plots are invaluable for con- 
densing and evaluating large amounts of data. These plots 
show, for example, that peak current is approximately propor- 
tional to wire length in the wire-length range from several 
meters to several hundred meters. Resistive loading of several 
ohms to several hundred ohms does not drastically affect peak 
currents or total energy. Calculations for reactive loads show 
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Fig. 30. Antenna pulse synthesis is straightforward using time-domain analysis [44]. Here we illustrate one possible procedure by its applica- 
tion to a conical-spiral antenna (the same one considered in Fig. 25). (a) The boresight radiation-field transfer function. (b) The transient 
voltage waveform required for the antenna to radiate the third derivative of a Gaussian pulse. (c) The actual radiated field produced by ap- 
plication of the voltage pulse in (a). 
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Fig. 31. (a) Wire-grid model of a light truck located over a perfect ground plane, used to compute the impedance characteristics of a 108-in 
whip antenna mounted as shown. (b) The antenna admittance demonstrates an input impedance of ~70 2 at 27 MHz, the frequency of 
operation for this emergency communication system (CB radio). 
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Fig. 32. EMP applications of time-domain solutions. A 10-m4ong wire 
of radius 6.738 X 10°? m and modeled with 30 equal-length seg- 
ments is illuminated from broadside by a nominal EMP pulse given by 
5.25 X 10* [exp (-4 X 10° t) - exp (~4.76 X 10° 2)] V/m. (a) The 
midpoint current and voltage that result for several midpoint load 
values. (b) Cumulative energy collected form the pulse by the wire, 
shown as a function of time (see (3a)). Observe that the collected 
energy is oscillatory in time, signifying the re-radiation of energy. 
(c) Energy absorbed by in the load as a function of time (see (3b)). 


that a single resistance in series with a resistance has little 
effect on total energy collected. 


IX. CONCLUSION 


In the above discussion, we have presented the rudiments of 
EM transient analysis via the use of direct time-domain integral- 
equation solutions, and demonstrated their utility for a variety 
of applications. We have emphasized the physical aspects of 
transient behavior and examined the numerical treatment in 
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Fig. 33. The value of a time-domain computation in finding an efficient 
early-time solution. (a) The peak current excited on a dipole by 
a broadside-incident EMP pulse (see Fig. 32) is plotted as a con- 
stant contour value as a function of wire length (m) and center load- 
resistance value (ohms). The time-domain computation permits the 
peak current, which occurs early in time, to be found without requir- 
ing a complete transient waveform. (b) Total load energy plotted as 
a constant contour value as a function of wire length (m) and center 
load-resistance value (ohms). 

The energy calculation requires that I7R decay to a small value be- 
fore terminating the computation. 


some depth, while introducing a minimum of mathematical 
detail. 

The single most important point we would like to have the 
reader retain is an appreciation for the practical utility of direct 
time-domain techniques in providing greater insight and under- 
standing of EM phenomena. Direct techniques also offer 
greater solution efficiency than transform techniques for many 
types of transient problems, the ability to handle nonlinearities, 
the convenience of wide-bandwidth information from a single 
calculation, the opportunity to use time-range gating to isolate 
interactions, and the possibility for directly obtaining the 
complex resonances of objects excited by EM sources. In 
closing, we hope the reader will be encouraged to employ direct 
time-domain computation and measurement where appropriate 
as an additional tool for solving EM problems. 


APPENDIX 


The derivation of (3c) follows readily from the approach used 
to develop the Poynting’s vector. From Maxwell’s equations 
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* Note that an explicit source (J) term is not included since the source 
here is an incident electric field. 
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and the H - product of (A1), the E- product of (A2), and their 
difference, we obtain 


V-EXH=-poH- E (A3) 


where we have used 

A'VXB-B:VXA=V‘AXB. 
If (A3) is integrated over a volume V that contains the con- 
ducting body, there is found 
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Upon using the Gauss divergence theorem and integrating with 
respect to time, we obtain 
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Now let the volume V be congruent with the volume (V,) of 
the conductor, assumed to be a wire, and extend a distance d 
beyond its surface. Then in the right side of (A5), we can 
approximate H and E by 
HIX A/a =X Al/2ar 
E © AQ/2n req 


(A6a) 
(A6b) 


since fH ~0 and $- E ~0 in the space between V and V,. 
Thus the right side of (A5) takes the form 
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where C is the wire contour and we have neglected the contri- 
butions to the volume integral of its ends (if the wire is open 
ended). When the r’ and ¢’ integrals of (A7) are performed, 
we finally obtain 
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which gives rise to (3c). 
The left side of (A5) can be similarly simplified to obtain 


t 
us ~ [ a { E-T dw’. 
-oo c 


Because E ‘J = Etan J = (Eexc + Escat)tan /=0 on a perfect 
conductor, equation (3a) for the collected energy immediately 
follows. When there is loss on the conductor, then Eygp/ = 
Ejossl = UR), from which (3b) also follows. 


(A9) 


PROCEEDINGS OF THE IEEE, VOL. 68, NO. 11, NOVEMBER 1980 


ACKNOWLEDGMENT 


The authors appreciate the assistance of discussions they had 
with several colleagues while preparing this work. These col- 
leagues include R. M. Bevensee, J. N. Brittingham, G. J. Burke, 
F. J. Deadrick, T. Lehman, and A. J. Poggio. They are also 
grateful to Margie Schmidt, Peggy Lorton, Debbie Koehler, 
and Sandy Johnston for their skill and diligence in typing the 
manuscript, and to Robert Waite for his patience and care in 
editing it. 


REFERENCES 


[1] E. K. Miller, “Some computational aspects of transient electro- 
magnetics,” Lawrence Livermore Lab., Livermore, CA, UCRL- 
§1276, 1972. 

A. J. Poggio and E. K. Miller, “Integral equation solutions of 
three-dimensional scattering problems,” in Computer Techniques 
for Electromagnetics, R. Mittra, Ed. New York: Pergamon Press, 
1973, ch. IV. 

G. Franceschetti and C. H. Pappas, “Pulsed antennas,” IEEE 
Trans, Antennas Propagat, vol. AP-22, p. 651, 1974. 

L. Felsen, Ed., Transient Electromagnetics. New York: Springer- 
Verlag, 1976. 

F. J. Friedlaender, Sound Pulses. 
Univ. Press, 1958. 

P. O. Brundell, ‘Transient electromagnetic waves around a cylin- 

drical transmitting antenna,” Erricsson Tech., vol. 16, no. 1, pp. 

137-162, 1960. 

O. Einarsson, “The step-voltage current response of an mfinite 

conducting cylinder,’? Trans. Roy. Inst. Technol. Stockholm, 

vol, 191, 1962. 

J, Rheinstein, ‘‘Backscatter from spheres: A short pulse view,” 

IEEE Trans, Antennas Propagat., vol. AP-16, p. 89, 1968. 

C. Flammer and H. E. Singhaus, “The interaction of electro- 

magnetic pulses with an infinitely long, conducting cylinder above 

a perfectly conducting ground,” AFWL Interaction Note, no. 144, 
1973, 

P. R. Barnes and D. B. Nelson, “Transient response of low fre-~ 
quency vertical antennas to high altitude nuclear electromagnetic 

pulse (EMP),” AFWL Interaction Note, no. 160, 1974. 

E. M. Kennaugh and D. L. Moffatt, “On the axial echo area of 
the cone sphere shape,” Proc. IRE (Corresp.), vol. 50, p. 199, 

1962. 

E. M. Kennaugh and D. L. Moffatt, “Transient and Impulse re- 

sponse approximations,” Proc. IEEE, vol. 53, p. 983, 1965. 

K. M. Mitzner, ‘‘Numerical solution for transient scattering from 

a hard surface of arbitrary shape—Retarded potential technique,” 

J, Acoust. Soc. Amer., vol. 42, p. 391, 1970. 

C. L. Bennett and W. L. Weeks, “Electromagnetic pulse response 

of cylindrical scatterers,’’ presented at G-AP Symp. Boston, MA, 

1968. 

See also, ‘“‘A technique for computing approximate electromag- 
netic impulse response of conducting bodies,’’ Purdue Univ., 

Lafayette, IN, Rep. TR-EE6811, 1968. 

E. P. Sayre, ‘Transient response of wire antennas and scatterers,” 
Elec, Eng. Dep., Syracuse Univ., Syracuse, NY, Tech. Rep. TR- 
69-4, 1969. See also “Time-domain radiation and scattering by 

thin wires,’’ Appl. Sci. Res., vol. 26, no. 6, pp. 413-444, 1972. 

C. L. Bennett and W. L. Weeks, “Transient scattering from con- 

ducting cylinders,” IEEE Trans. Antennas Propagat., vol. AP-18, 

pp. 627-633, 1970. 

E. K. Miller, A. J. Poggio, and G. J. Burke, “An integro-differen- 

tial equation technique for the time-domain analysis of thin-wire 

structure. Part I, the numerical method,” J. Comput. Phys., 

vol. 12, no. 1, pp. 24-48, 1972; “Part II, Numerical results,” J. 

Comput. Phys., vol. 12, no. 2, pp. 210-233, 1972. 

A. J. Poggio, “Space-time and space-frequency domain integral 

equations,” IEEE Trans. Antennas Propagat., vol. AP-19, pp. 

702-704, 1971. 

T. K. Lui and K. K. Mei, “‘A time-domain integral equation solu- 

tion for linear antennas and scatterers,” Radio Sd., vol. 8, nos. 

8-9, pp. 797-804, 1973. 

T. T. Wu, “Transient response of a dipole antenna,’’ J. Math. 

Phys., vol. 2, no. 6, pp. 892-894, 1961. 

L. B. Felsen, ‘‘Diffraction of the pulsed field from an arbitrarily 

oriented electric or magnetic dipole by a perfectly conducting 

wedge,” SIAM J. Appl. Math., vol. 26, no. 2, p. 306, 1973. 

L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves. 
Englewood Cliffs, NJ: Prentice-Hall, 1973. 

V. H. Weston, “Pulse return from a sphere,” IRE Trans. Antennas 


[2] 


[3] 
[4] 
[5] 
[6] 


London, England: Cambridge 


(7] 


[8] 


[9] 


[10] 


(11] 


[12] 
(13) 


[14] 


(15] 


[16] 


[17] 


[18] 


[19] 


[20] 
(21] 


[22] 


[23] 


MILLER et al.: TRANSIENT RADIATION AND SCATTERING FROM WIRES 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


Propagat. (Special Supplement), vol. AP-7, pp. S43-S51, Dec. 
1959, 

D. L. Sengupta and Y-P Liu, ‘‘Analytical investigation of wave- 
forms radiation by a resistively loaded linear antenna excited by a 
Gaussian pulse,’ Radio Sci., vol. 9, no. 6, pp. 621-630, June 
1974, 

C. E. Baum, “Emerging technology for transient and broad-band 
analysis and synthesis of antennas and scatterers,’’ Proc. IEEE, 
vol. 64, pp. 1598-1616, Nov. 1976. 

L. Marin, ‘“‘Natural-mode representation of transient scattering 
from rotationally symmetric bodies,” JEEE Trans. Antennas 
Propagat., vol. AP-22, p. 266, 1974. 

K. R. Umashankar, T. H. Shumpert, and D. R. Wilton, “‘Scatter- 
ing by a thin-wire parallel to a ground plane using the singularity 
expansion method,” JEEE Trans. Antennas Propagat., vol. AP-23, 
p. 178, 1975. 

M. L. Van Blaricum and R. Mittra, ‘A technique for extracting 
the poles and residues of a system directly from its transient re- 
sponse,”’ [EEE Trans. Antennas Propagat., vol. AP-23, p. 777, 
1975. 

E. K. Miller, F. J. Deadrick, H.G. Hudson, A. J. Poggio, and J. A. 
Landt, ‘“‘Radar target classification using temporal mode analysis,” 
Lawrence Livermore Lab., Livermore, CA, UCRL-51825, 1975. 
J.N. Brittingham, E. K. Miller, and J. L. Willows, ‘‘Pole extraction 
from real-frequency information,’ Proc. IEEE, vol. 68, pp. 
263-273, Feb. 1980. 

E. K. Miller, J. N. Brittingham, and J. L. Willows, “Identification 
of E. M. spectrum by known pole sets,” Electron. Lett., vol. 13, 
no. 25, p. 774, 1977. 

C. L. Bennett, A. M. Auckenthaler, R. S. Smith, and J. D. DeLo- 
renzo, ‘‘Space-time integral equation approach to the large body 
scattering problem, Rome Air Development Center, Rome, NY. 
Rep. RADC-TR-73-70, 1973. 

A. M. Nicolson, C. L. Bennett, D. Lamensdorf, and L. Susman, 
“Applications of time domain metrology to the automation of 
broad band measurements,’’ JEEE Trans. Microwave Theory 
Tech., vol. MTT 20, pp. 3-9, Jan. 1972. 

R. M. Bevensee, F. J. Deadrick, E. K. Miller, and J. T. Okada,” 


[35] 
[36] 
[37] 


[38] 


[39] 
[40] 


[41] 


[42] 
[43] 


[44] 


[45] 


146] 


[47] 


1423 


Validation and calibration of the LLL transient-electromagnetic 
measurement facility, Lawrence Livermore Lab., Livermore, CA, 
UCRL-52225, 1977. 

F. Rohrlick, Classical Charged Particles. 
Wesley, 1965. 

W. K. H. Panofsky and M. Phillips, Classical Electricity and Mag- 
netism. Reading MA: Addison-Wesley, 1956. 

J.A. Landt, E. K. Miller, and ML. Van Blaricum, ‘“WT-MBA.LL1B: 
A computer program for the time-domain response of thin-wire 
structures,’’ Lawrence Livermore Lab., Livermore, CA, UCRL- 
51585, 1974. 

E. K. Miller, F. J. Deadrick, and J, A. Landt, ‘“Time-domain anal- 
ysis of non-linear loads,’’ Lawrence Livermore Lab., Livermore, 
CA, UCRL-52172, 1976. 

R. F. Harrington, Field Computation by Moment Methods. New 
York: Macmillan, 1968. 

R. Mittra, Ed., Numerical and Asymptotic Techniques in Electro- 
magnetics. New York: Springer-Verlag, 1975. 

D. E. Merewether, “Transient currents induced on a metallic 
body of revolution by an electromagnetic pulse,” JEEE Trans. 
Electromagn. Compat., vol. EMC-13, pp. 41-44, 1971. 

C. L. Bennett and G. F. Ross, ‘‘Time-domain electromagnetics 
and its applications,” Proc. IEEE, vol. 66, no. 3, Mar. 1978. 

H. K. Schuman, “Time-domain scattering from a non-linearly 
loaded wire,” IEEE Trans. Antennas Propagat., vol. AP-22, 
p. 611, 1974. 

T. K, Liu and F. M. Tesche, “Analysis of antennas and scatterers 
with non-linear loads,’’ IEEE Trans. Antennas Propagat., vol. AP- 
24, p. 131, 1976. 

T. K. Sarkar and D. D. Weiner, “Scattering analysis of non-linearly 
loaded antennas,’’ JEEE Trans. Antennas Propagat., vol. AP-24, 
p. 125, 1976. 

J. A. Landt and E. K. Miller, ‘‘Transient response of the infinite 
cylindrical antenna and scatterer,” JEEE Trans. Antennas Propa- 
gat., vol. AP-24, p. 246, 1976. 

E, K, Miller and J. A. Landt, ‘Short pulse characteristics of the 
conical spiral antenna,’’ JEEE Trans. Antennas Propagat., vol. 
AP-25, pp. 621-630, Sept. 1977. 


Reading, MA: Addison- 


Fig. 14. Physical and numerical modeling errors can be assessed using 
models like these. Shown here are a scale-model and pipe and thin- 
wire models of a Boeing 747 aircraft [34]. 


